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OPTIMAL MODAL CONTROL OE POWER SYSTEMS 

Recent decades have witnessed a phenomenal growth 
in the size and interconnection of power systems all over 
the world. The problems of design and operation of these 
systems are complex and require sophisticated techniques, 
especially from the area of system theory, for their 
solution. 

Production of large amounts of power at low cost 
has necessitated employment of generators of high capacity, 
exploitation of remote hydro-electric sources and trans- 
mission of bulk power by e.h.v, transmission lines, and 
increased interconnections to obtain economy and reliabi- 
lity. The fundamental problems in the design and operation 
of such systems are those associated with system stability, 
economic operation, reliability of supply to the consumers., 
security and operation within the constraints such as those 
on system voltages and frequency. Many of these problems 
can be posed as problems in the area of automatic control. 
Consequently there has been considerable activity, in 



recent years, in formulating these problems and determining 
efficient solution techniques to them using techniques of 
modem control theory. 

Earlier approaches to improve the power system 
stability by the employment of fast acting exciter voltage 
regulating 'systems resulted in improved transient stability 
of the system. However, these systems generally exhibited 
poor damping of machine angle swings. These swings cause 
variation in system voltages and frequency and in extreme 
cases may even result in system instability. An appropriately 
designed feedback controller may overcame this problem. 
Therefore, in the recent past, researchers have paid great 
attention towards the development of suitable controls to 
damp out quickly the machine angle oscillations and at the same 
time retain the improvement in the transient stability. A 
major difficulty with this problem is that the system 
equations are nonlinear and of high order. An approach 
which is usually employed makes use of a linear system 
model valid for small departures of the system from its 
nominal operating conditions. Optimal infinite time linear 
regulator theory is then used to obtain linear , time- 
invariant feedback laws. The following difficulties are 
encountered in this approach* (a) It is difficult to 
determine apriori what the state and control weighting 
matrices, Q and R respectively, should be to obtain 



satisfactory closed loop system performance. It is 
necessary to solve the optimal control problems successively 
fo** various Q and R matrices, till satisfactory system 
performance is obtained. The method does not provide the 
designer with a 'feel’ for the problem, and is computa- 
tionally unattractive, (b) The method requires accessibi- 
lity of all states of the system for measurement. The use 
of observers which overcomes this drawback, introduces its 
own computational complexity. (c) The system should be 
completely state controllable. (d) The feedback law 
obtained for one operating point of the system may be 
inappropriate at others. (e) Furthermore, there is a 
problem associated with feedback gains. From the practical 
point of view, large feedback gains should be avoided. The 
approach referred to above does not provide direct control 
over the feedback gains . 

This thesis attempts to overcome some of the problems 
mentioned above by developing suitable techniques stemming 
from modal control theory. Like the classical root locus 
technique, modal control approach provides the designer with 
a ’feel* for the problem, and his engineering judgement can 
be better utilized for developing suitable designs. In many 
practical situations, system behaviour is governed mainly by 
a small number of 'dominant' eigenvalues. A suitable design 
can usually be obtained by determining a feedback controller 



which shifts these eigenvalues far enough into the left 
half plane without disturbing the locations of the non- 
dominant eigenvalues. Complete state controllability is 
not necessa.ry for achieving this objective. This thesis 
reports a technique for minimizing feedback gains to 
obtain optimal modal controllers and procedures which 
facilitate determination of modal controllers, optimal 
and non-optimal, for high order systems. In the case of 
high order systems, grouping of eigenvalues and use of 
appropriate control inputs to shift these eigenvalues is 
proposed. Using these techniques and coupled with the 
knowledge regarding sensitivity of the eigenvalues to 
system parameter perturbations, wide range modal control- 
lei's have been developed. These controllers give 
satisfactory system performance for a range of operating 
conditions. The thesis also considers design of incomplete 
state feedback and output feedback modal controllers. Modal 
control technique developed has been also applied to design 
an output feedback dynamic controller. 

The design techniques developed pertain to finite- 
dimensional linear, time -invariant systems; they can 
therefore be applied in situations where the physical 
systems can be adequately represented by such mathematical 
models. In addition to applying the technique to the 
control of a power system with its associated controls, the 



thesis also illustrates the use of some of the techniques 
to the design of Load Frequency Controllers and controllers 
for high order systems. Specifically modal controller 
design of a forty-first order system describing a chemical 
plant has been given for the purpose of illustration. 

A chapterwise summary of the work reported in this 
thesis is given below: 

The problem of design of feedback controllers for 
power system is introduced in Chapter 1, A review of the 
recent developments in this area is given. 

Chapter 2 is devoted to the development of the 
state space model for a power system which is used for 
illustrating the techniques of control developed in 
Chapters 3» 4, 5 and 6. The power system considered 
consists of a synchronous generator connected to an 
infinite bus through a transmission line and has an 
exciter-voltage regulator and a turbine speed governor. 

The generator is -driven either by a hydraulic turbine, or 
a steam turbine (as a prime-mover). The nonlinear 
equations describing the system dynamics are derived. 
Linear models are obtained by linearizing the nonlinear 
equations. The model for the design of Load Frequency 
Control scheme for a two area system is also briefly 
described. The control problems associated with the 
above models have been enunciated. 



Chapter 3 introduces the concept of optimal modal 
control by suggesting appropriate performance indices. The 
indices involve minimization of feedback gains. The 
closed loop eigenvalues are decided as in classical control 
theory based on ’feel’ of the problem. A solution technique 
for designing an optimal modal controller is then described. 
Using this algorithm optimal modal controllers have been 
designed for a power system. The algorithm has also been 
applied to the design of a load frequency controller for 
a two area problem. The responses are compared with con- 
ventional schemes as well as uncontrolled responses. 

The problem of designing controllers for high order 
systems are discussed in Chapter 4. A straightforward 
application of the technique developed in Chapter 3 will 
result in high feedback gains. This will be so in spite 
of the feedback gain optimization, due to the constraints 
imposed on the structure of the feedback matrix. By 
relaxing the constraint, an alternative design procedure 
for optimal modal control of large systems has been 
developed. This procedure is based on grouping of eigen- 
values and moving these groups to desired locations 
successively through appropriately chosen inputs. Modifi- 
cation of this procedure which results in advantages in 
computation as well as implementation has been indicated. 
Three illustrative problems for the design of controllers 
for power systems have been solved .and the results discussed. 



To demonstrate the efficacy of the algorithms developed, 
they have been applied to the controller design for a 
chemical plant having fortyone state variables and eight 
control inputs. 

Linear regulators for power systems are usually 
designed fo** a particular set of system operating condi- 
tions. However, these conditions change with the load 
demand on the system. In Chapter 5 , the design of a wide- 
range modal controller is developed using modal control 
theory and eigenvalue sensitivity analysis. The design 
reduces the variation in closed loop pole locations with 
change in operating conditions. A numerical example of 
the design of the wide range controller has been given and 
the practical implementation discussed. 

Modal controller design requires feedback of all 
the state variables. In practice it may be difficult or 
impossible to do so. A technique to design modal con- 
troller with accessible state or output feedback is 
described in Chapter 6. The chapter also includes the 
design, using modal control approach of a dynamic output 
feedback controller. 

Chapter 7 reviews the work reported and contribution 
made in this thesis and concludes with an assessment of the 
scope' for further work in this area. 



CHAPTER 1 


INTRODUCTION 

1*1 POWER SYSTEM CONTROL PROBLEM 

The phenomenal growth in the size of power systems 
and the increase in interconnections to improve reliability 
of bulk power supply results in system stability (both 
transient and dynamic) problems* By transient stability 
we imply system response to large disturbances and dynamic 
stability pertains to small signal behaviour around an 
operating point. Both phenomena are of concern to system 
planning and operating people. Construction of super 
thermal stations, commissioning of large capacity genera- 
ting units, use of extra-high- volt ages for large distance 
transmission, hydraulic generation at remote locations and 
the design of modem alternators with large synchronous 
reactances reduce the stability margin. To improve the 
transient stability of the power system, various control 
strategies have been investigated [l-8] . Considerable 
improvement in the transient stability has been obtained 
by using fast acting excitation systems [2]. Though the 
damping of rotor angle oscillations is initially more, it 
becomes poor after the first cycle. It is also known that 
modern fast acting AVR's and governors while contributing 
to improved steady state response and dynamic control can 
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aggravate the problem of small signal dynamic stability. 
The current practice of running the gerte^ators at high 
power-factor and the generation of reactive power due to 
adoption of extra-high voltage transmission lines and 
cables, requires that more emphasis be placed on the deve- 
lopment of compensating controls to stabilize the system. 

Power system simulation studies are used to extend 
the simulation of power system behaviour for several 

minutes after the first swing and provides the transition 

- 

between the first swing transient stability which gives 
immediate post-fault behaviour and the load flow which 
describes the steady state of the post-fault conditions. 

A properly planned and operated power system permits 
satisfactory transient performance following a disturbance 
and then returns to a desirable steady state. Stability 
studies are done to see that the system meets the require- 
ments [9,10]. 

After a disturbance oqcuts on the system, the state 
of the system such as voltage levels, frequency, tie line 
power flows and line currents will change. In extreme 
cases, the synchronous machines in the system may become 
unstable. Control of a power system to bring back the 
system state to normal operating condition is an essential 
requirement [ll]. Moreover, the control actions must be 
such as to restore the state to the normal condition in 
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the shortest possible time after a disturbance occurs on 
the system. 

In the methods of power system control the deviation 
in the states such as voltage, frequency and machine angles 
are detected by means of sensors. Depending on the error, 
the voltage regulator and excitation system and the 
governor and turbine system take corrective steps to keep 
the system states within limits. In recent year's, the 
effect of excitation and speed governing systems on power 
system performance has been studied extensively. Improve- 
ments in the system performance through subsidiary feed- 
back applied to the basic voltage or frequency regulating 
system have been investigated. The feedback of supple- 
mentary signals obtained from rotor angle, rotor speed and 
terminal voltage into the excitation system has been shown 
to be very effective [3»12,13]» Classical control tools 
such as Routh-Hurwitz, Nyquist and Root locus technique 
have also been used to analyse the dynamic stability 
phenomena of a power system [14-16]. However, these 
techniques are cumbersome to apply even to a simple power 
system and almost impossible for large multimachine inter- 
connected power systems. 

1.2 MODERN TRENDS AND THE STATE OP ART 

The availability of digital computers has made 
it possible to undertake stability studies of large 
interconnected power systems. Considerable interest 
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has been shown in recent years [6-8,17-47] in the applica- 
tion oi the restilts of modem control theory to the 
optimization of power system performance. The application 
of the modern control theory has enabled the optimization 
to be carried out in a systematic and logical manner. The 
philosophy of system optimization requires that the per- 
formance index (J ) in the form of the integral of a 
function of the system variables, x> and the control 
functions, u, be defined and minimized over a specific 
time interval. The state variable transients and control 
actions are minimized simultaneously. 

The dynamics of the power system are represented 
by a set of nonlinear differential equations 

x = f (x, u, t) (1.1) 

The optimal control law is obtained by solving a set of 
differential equations obtained by applying Pontryagin’s 
minimum principle. Some researchers [17-19,48] attempted 
to solve the nonlinear power system optimization directly. 
In references [17] and [18] the optimal control has been 
obtained as a function of time, using a N ewt on-Raphs on 
gradient minimization and quasilinearization technique 
respectively. Rajagopalan and Hariharan [48] proposed 
a bang-bang form of excitation input and the switching 
time is calculated to minimize a performance index by a 
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hill climbing method. The open loop form of control 
calculated by the above procedures As- highly sensitive 
to variation of system operating conditions. Quasi-Fewt on 
method of function minimization to derive closed loop 
constant gain feedback control to optimize nonlinear 
system response has been attempted in reference [l9]« 

In this approach, the nonlinear equations (l.l) are to 
be solved to evaluate the performance index. 

The solution of nonlinear eqns. (1.1) is computa- 
tionally difficult. However, for studying the transients 
due to small disturbances, the original nonlinear system 
can be linearized about the operating point to the linear 
matrix form 

* 

x = A x + B u (1.2) 

where the variables now denote their deviation from 
their nominal operating conditions* 

A linear regulator with constant feedback gains, 
can be designed, using modern control theory, to obtain 
satisfactory closed loop system performance. There are 
mainly two approaches to the design of linear regulators: 
(1) procedures based on linear optimal regulator theory 
and (2) procedures based on pole assignment techniques. 

1.2.1 Optimal Control of Synchronous Machine: 

Most researchers [6-8,20-36] have obtained the 
optimal control for a power system through the application 
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of linear regulator theory which involves the solution of 
the matrix Riccati equation. The optimal control law so 
obtained requires the feedback of all the state variables. 
The direct measurement of all the state variables may not 
be possible in a system. Even though an observer may be 
constructed to estimate the states, it may not be possible 
to implement it. In a realistic power system, to provide 
feedback of certain measurable states, may not be practi- 
cable. For example, in a hydro system the equipment to 
measure the change in water level is situated far away from 
the generating station. Transmission of such signals may 
be expensive and noise is introduced in transmitting such 
signals. These constraints prompted the design of 
controllers based on the measurement of available states 
or outputs [ 21 - 25 ] • Davison and Rau [ 21 ] obtained feed- 
back laws in terms of the various combinations of outputs 
using parameter optimization method and the solution of 
Liapunov’s matrix equation. DeSarkar et al [ 22 ] and 
Pai et al [39] designed output feedback controllers by 
detecting the feedback from undesirable states. However, 
controller designs based upon above approach may not give 
satisfactory response and in extreme cases the system may 
even become unstable. In reference [ 23 ] a suboptimal 
control policy obtained through a reduced model which 
contains only dominant eigenvalues of the original system 
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is proposed. Habibullah et al [24] obtained the output 
feedback controller by choosing a state space model having 
outputs as state variables. Eaina et al [25] designed a 
feedback control system based on quantities easily measur- 
able at the synchronous generator locations. Furthermore , 
the designs of references [24*25] have been shown to give 
satisfactory response for a wide range of operating condi- 
tions. 

The controllers designed by linearizing the system 
equations about a particular operating point may not give 
satisfactory response for other operating conditions. An 
additional correcting feedback signal to compensate for 
small variation in operating conditions has been suggested 
[ 26j * Optimal regulator theory incorporating the trajec- 
tory sensitivity function in the performance index has been 
suggested by Elmetwally and Rao [27*50,32] to design low 
sensitivity excitation control. Design based on linear 
optimal regulator theory and eigenvalue search technique 
has been also suggested [24], to design wide power range 
controllers for power systems. 

1.2.2 Load Frequency Control Problem; 

Ever since the pioneering work of Elgerd and 
Fosha . [37»38 ]» there has been continuing interest in 
the multiarea load frequency control problem [39-45] • 

Venkat et [.41] designed load frequency controller 
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for systems with -unknown disturbances utilizing the state 
augmentation approach of Smith and Davison*. Ramamoorty 

et al [43] and Rao et al [44] designed regulators for two- 
area power systems via Liapunov's second method and 
utilizing minimum settling time theory. The importance 
of the dominant time constant of the closed loop system 
in designing the regulators has been emphasized by them. 
Furthermore, Rao et al [44] suggested a criterion for 
choosing the frequency bias setting. In a recent paper 
Elmetwally and Rao [45] designed a controller with assumed 
constrained structure to study the decentralized control 
of multiarea load frequency control problem. Their design 
of local controller is easily impl ementable and gives 
comparable results with those obtained by centralized 
controller. 

1.2.3 Pole Assignment Procedures: 

Recently in the literature, emphasis has been 
given to the locations of the closed loop system poles, 
while designing the regulators for power system control 
[24,31»36,49-5l] * Elangovan et al [ 36 ] suggested a method 
to design controllers with prescribed eigenvalues. In 
their approach, the state weighting matrix in the perfor- 
mance index is selected to give desired closed loop eigen- 
values locations. Using pole assignment procedure, Seraji 
[ 31 ] designed controller for a single input synchronous 
generator. Pai et al [49] applied modal control theory 
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for designing controllers for a power system with two 
inputs. However, their technique utilizes only one of 
the two available inputs to shift one or more eigenvalues 
at a time. 

1.3 OBJECTIVE AND SCOPE OE THE THESIS 

1. In the optimal linear regulator theory, it is 
necessary to systematically change the state and control 
weighting matrices, Q and R respectively, in the perfor- 
mance index ( J ) , till a satisfactory closed loop system 
behaviour is obtained. It is therefore better to employ 
design techniques which directly yield :-.l.thb.. required 
closed loop system eigenvalue locations. This would avoid 
the problem of having to iteratively solve a sequence of 
optimal linear regulator problems. Furthermore, the 
optimal linear regulator theory does not provide the 
designer with a *feel f for the problem, and is computa- 
tionally unattractive, especially for high order systems. 
Pole-placement techniques coupled with a performance 
criteria for optimization is a desirable goal to achieve. 

2. Controller designed for power systems using the 
optimal infinite time linear regulator theory or modal 
control techniques [52-55] requires feedback of all the 
states of the system. It is seldom that all of them can 
be monitored directly in the case of power systems. The 
use of observers ^though overcomes this problem, introduces 
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its own computational complexity. Furthermore, the require- 
ments that the system should he completely controllable for 
designing optimal controller and completely observable for 
the design of observer are restrictive. In certain systems 
the output variables which are linear combination of state 
variables may only be available for measurements. Hence 
controllers based on output feedback are practically 
desirable. 

3. In the case of large scale power systems, there 
is a problem associated with the feedback gains. The 
complexity of the problem increases in arranging feedback 
to all the control inputs. Excessive computational efforts 
may be required in optimizing the controller design for 
large systems with more number of control inputs. A 
technique in moving eigenvalues in groups using one or 
more selected inputs which at the same time minimizes 
feedback gains will result in an improved controller. 

4. The linear regulators (optimal and modal 
controllers) are designed for a particular operating 
conditions. In a power system, the operating conditions 
do change. Therefore, the controllers designed for one 
operating condition may be inappropriate for other 
operating points. The need for a wide range controller 
is therefore clear. 
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Tiie thesis addresses itself to the problems enumerated 
above. Accordingly the ob jectives .of the thesis are set forth 
as follows . 1 

(i) to develop suitable techniques stemming from 
modal control theory for the design of modal controllers 
for power systems so as to minimize the feedback gains. The 
techniques should be such that the designer’s engineering 
judgement can be better utilized for developing suitable 
designs. 

(ii) to suggest a suitable procedure using eigen- 
value grouping technique for minimizing feedback gains to 
obtain optimal modal controllers. Computationally simple 
procedures for the design of modal controllers, optimal and 
non-optimal for high order systems are to be developed. The 
possibilities of designing controllers with feedback to 
small number of inputs without appreciably affecting the 
feedback gains are explored. 

(iii) to propose control laws using sensitivity 
techniques so that, the effect of the wide range change in 
operating conditions is small on the locations of the closed 
loop pole locations. 

(iv) to present design techniques for designing 
physically realizable modal controllers with easily acces- 
sible states or output feedback. 

The chapterwise summary of the various chapters in 
this thesis are given in the synopsis. 



CHAPTER 2 


MODELLING OE POWER SYSTEMS 

2.1 INTRODUCTION 

An adequate mathematical model is a prerequisite 
for any application of modern control theory to design 
controllers for power systems. With the availability of 
fast and large digital computers, sufficiently complex 
power system models can be analysed. In general power 
system model is nonlinear, comprising of the differential 
equations of the synchronous machine and its controls and 
the algebraic equations governing the transmission network. 
Application of optimal control theory to such models is 
difficult. Therefore the usual practice is to linearize 
the model around an operating point and then develop the 
control laws. Linear models are in wide use due to the 
ease with which one can apply the results of multivariable 
control theory. High order models yield better accuracy 
but difficulties arise in computation and practical imple- 
mentation, of control strategies. It is therefore necessary 
to strike a balance between ease of implementation and 
engineering accuracy which is desirable. The literature 
in the area of dynamic modelling of power systems is quite 
profuse [56-65]. Hence it is not proposed to discuss the 
modelling aspect in detail. 
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The basic power system components to be modelled 
are (i' transmission line, (ii) synchronous machine, 

(iii) voltage regulator and excitation system and 

(iv) primemover and governor system. In dynamical studies 
the variation in the rotor angles is essentially a low 
frequency phenomenon. Hence the transients in the trans- 
mission network are neglected and steady state phasor 
representation of the transmission network is used. The 
synchronous machine is represented by Park's equations 
[56]. To reduce the complexity of the model, two simpli- 
fications are made (i) the stator transients are neglected 
and (ii) the speed voltage terms are evaluated with w equal 
to w Q (synchronous speed). The first assumption is consis- 
tent with the representation of the transmission network. 

It is further assumed that the resistance of the trans- 
mission line and the effect of magnetic saturation in the 
machine are neglected. The voltage regulator and the 
exciter system is represented by a single time constant. 

The primemover and governor system is represented by two 
time constants [65], which is simplified to a single time 
constant in the case of fast acting governor and pilot 
valve system. 

Various aspects of modal control theory developed 
in this thesis have been applied to basically three dyna- 
mical models of the power system. (1) Power system 
considered by Yu et al [6], (ii) Power system considered 
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by Davison and Rau C 21 ] and (iii) Two area load frequency 
control system considered by Fosha and Elgerd [ 38 ] . 

A brief description, the performance equations, 
the linearized model and the data for the above systems 
are given in the subsequent sections. The details are 
contained in the appropriate references. 

2.2 POWER SYSTEM MODELS 

2.2.1 Power System Considered by Yu et al [6]: 

The power system considered consists of a synchronous 
generator driven by a hydraulic turbine, having an exciter- 
voltage regulator and a turbine speed governor. The genera- 
tor is connected to an infinite bus via a long transmission 
line [Figure 2.1]. 

The increment in field voltage i.e. A v f is 

controlled through a signal u^. This signal is applied 

througa the voltage control feedback loop having transfer 

function [p. /(l + t p)], to the summing junction of the 
s s 

exciter-voltage regulator system [Figure 2.2]. 

Increment in the power input AP^ to the machine 
is controlled by a control signal u^* This control signal 
is applied via a servomotor represented by a simple transfer 
function [u /(I + t p)] to the summing junction of the 

ci cl 

governor system [Figure 2.3]. 

The synchronous machine has a three phase balanced 
winding and the resistance of the winding is smaller than 




FIGURE 2.3 GOVERN OR-HYDRAU1I G SYSTEM 
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its reactance. The rotor carries the field winding only. 
The machine terminal conditions (Figure 2.1) give in 
d-q components the relation, [6] 


(v, + jv ) (jB ) + [(v, + jv ) - (v Sin6 + jv Cos6)]/jx 


\ ^ .- 0 — — . «- 0 - 

* % + ^q 


( 2 . 1 ) 


The park's equations to he considered are 


v 

“fq 

% 

where 



“*d 


-x i 


q q 


1+ ^o p 


1+ ^o P 



( 2 - 2 ) 

(2-3) 

(2.4) 

(2.5) 


v Q = infinite bus voltage 
v^ = a field voltage 

v d ,v^ = d and q component of v^, respectively 
i ,i = d and q component of armature current i, 

'"** V£ 

r espe ct ively ? , 

= d-axis and q-axis flux linkages k respectively 

^ A 

6 = torque angle (radians) 
u> 0 = synchronous speed (rad/sec.) 
u> - speed (rad /sec. ) 

x^,x d , t^ o , , x^ = appropriate synchronous machine 


constants 
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let 

x *= p t (2.6) 

where p is the time scaling factor. Then 

& 


n4 

Aw 

P 

(2.7) 

d A 6 
dT 

n 

(2.8) 

dn = 
dT 

E A P ± - p D n - A p ] 

P 2 M 1 e 

(2.9) 


where 

AP^ = increment in the mechanical power input 
AP g = increment in the power of the electro- 
mechanical energy conversion 
D = damping coefficient. 

The increment in the mechanical power input is 
approximated by [6], 


AP i = A g + 1.5 A h (2.10) 

where Ag and Ah are the increments in the gate opening 
and the water head respectively-. The power of electro- 
mechanical energy conversion is given by the expression 


v ¥ Sin 6 

P = -2— £ + 

e 

d do 


(x* - x ) 

— ~ SL- v* Sin 26 

2 Y» 1 0 

^ X d q 


where 


?e " P o 


( 2 . 11 ) 


P Q = the power delivered by the machine at the 
operating conditions. 
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linearizing eqn.(2.11) around the o 
v Cos 6 (x* - x ) 

ap « [ ~° - ^ g a 4- -A— _a_ 


operating point, one gets 


Z d *Jo 

V Sin 6 
+ -2 : AY- 

z d T ao 


-A- v Cos 26 ] A 6 

X' Z 0 °- 

d q 


(2.12) 


Substituting eqns. (2.10) and (2.12) in eqn.(2.9)> one gets 

ft - -f t [- . iifcil T . Oos 2 3 4 6 

Cl T o 2m -tT » _ « y | yp O O 

X d A q 


x a T ao 


v Sin 6^ 

-pin- -2 2 AY + Ag + 1.5 Ah} (2.13) 

X’ t‘ 1 

d do 

The rate of change of flux linkages of the field winding 


is given by 

d Y v X.'Ff (x - x») 

— £ = —X - -X - Jl— v Cos 6 

P P*# 5 *, P x a ° 


(2.14) 


On linearization around the operating point, the above 
equation yields 


2J2e. - V J-- - ^ 1 31 V Sin 6 46 

4T P P z a T ao P H ° (2 . 15) 


where 


x! = direct axis transient reactance of the 
d 

synchronous machine 
X* = (1 - x B) x£ + x 
X„ ^ (1 - x B) x 4- x 

q q 
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A , 

Z d = (1 - x B) x d + x 

B = transmission line susceptance (Figure 2.1). 


The exciter voltage regulator system equations in 
state variable forms are obtained from Figure 2.2 as 


dA v f 
dx 

= - -ill + 

jrf e <‘ 4 T t- 4 V 

(2.16) 

dA v g 

_ . f • . . 

is_„ 

(2.17) 

dv 



where 




< 

Cvj CT 1 
> 

+ 

ii 



and the 

first order variation in v^. is 

approximated as 

Av t 

- ^2 4 v. ♦ 

T to 4 

V 

T to 1 

(2.18) 


On substituting for Av, and Av , the above equation 
reduces to 


Av t = 


v xv, Cos 6 x' v Sin 6 
0 .[ o _ d qo o j A 


x v 

+ S£ 


*<! 


v xjp- 

to d do 


(2.18a) 


Substituting (2.18a) in (2.16), the eqn.(2.16) reduces to 


d A v f 

d*c 


A v 

FT 


f 


pT 


v x v,„ Cos 6^ x ' v 

{ 2_ [_a d0 0 d 


e 


'to 


■v 




SinS 


x; 


X V 


3SL 


to *d T d0 


A% 


- A v 


s} 


(2.19) 


•]A6 
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The primemover governor system equation in state 
variable form are obtained from [Figure 2.3] 


d A g 
dT 

O « 1 

' fi'e ® + pT g 


Agf 

“a , 

dT 

pT a + 

P T a 

dA h 
dT ~ 

2d Ag 

dT 

2_ 

Pt 
r w 


( 


-im- id 


w o 


u^ 


Ah 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


The state variables selected for the linearized 
model are 

iT 


x.. = [fi 6, n,A¥ f » Av f , Av g , A g, Ag f , Ah]' 


(2.23) 


and the linearized power system model is given by the 
eqns. (2.8), (2.13), (2.15), (2.19), (2.17) and (2.20) to 
(2.22) which are put in standard state space form as 

S = A x + B u (2.24) 

where 


£ = [u v , u g ] T ( 2 . 25 ) 

Li ^ 

(2 - 26) 

K s 
K u o 

V -ff (2 - 27) 

St 

V = voltage control feedback loop gain 
0 

p = governor actuator gain 

8L 
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x — voltage control feedback loop time constant 
t = governor actuator time constant 

£l 


A = 


~o 

1 

0 0 

0 

0 

0 

0 


-D 

q n 

0 

1 

n 

l.p 

22 

PM 

23 

p 2 M 

V/ 

p 2 M 

a 31 

0 

1 

a 33 F 

0 

0 

0 

0 


0 

-1 

~»e 

0 

0 

0 

a 41 

*« Px e 

^ x e 

0 

0 

0 0 

-1 

0 

0 

0 

0 

-1 

A 

A 

-a 

-1 

A 

“o"g 

u u 

u 



u 

0 

0 

0 0 

0 

0 

-1 

? T a 

0 

r 

o 

2 

0 0 

A 

20 

2 

-2 

0) _T 

0 g 

U 


P T 

g 

P W 

(2.28; 

re 

• 

~ 1 ' * [ 

v o Cos 6 0 

W 

f0 + 

( M * 

to 

o 

o 

cm o 

i> 

2 u 

a 2l = 

p 2 M 

H Mo 


X 'd X , 


v Sin 6 
o o 





a 23 = 

p 2 M 

*4 Mo 






- ii. 

- x') 

d TT 

Sin 6^ 
0 




a 31 = 

P 

x* 0 

*d 





X, 


a 33 " ~ 


* X d T do 



22 . 


a 41 = 


V r, x v, ^ Cos 6^ 
_ e o j- q do o 


p-r_ v 


e ’to 


X. 


x a v a° an 6 0 - 


H 


a 43 " 


M = 


T - 


e 


x v 


qo 


Xj t ! 


P T e v to ~d v do 


inertia constant 
exciter time constant 


(2.29) 


t = gate time constant 
S 

t = water time constant 
w 

|i = exciter gain 
© 

o = permanent droop. 


The control matrix 

”o 0 0 0 


B = 


1 

0 0 


0 

0 


0 

1 


0 

0 


“I T 


(2.30) 


The relevant data for the power system considered 
above are, [6] ; 


Transmission Line 

x = 0.7417, B = 0.1339 

Synchronous Machine .. : 

x, = 1.00, x* = 0.27, x„ = 0.60, t ' = 9.00, 

d d q do 

M = 0.2122, D = 0.00537. 

Exciter and Voltage Regulator System 

= 10.00, = 1.0, -c g = 0.50. 
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Governor and Hydraulic System 

o = 0. 045 > Tg = 0 . 10 , = 0.01, - 1.60* 


The numerical values of u 

s 

are included in u and u , 

v g 


and u are not needed as they 
a 

respectively. 


Nominal Operating Conditions ; 

( ■■>, •' 1 

= 0.735 p »u. , Q o = 0.034 p.u. , v^ 0 = 1.05 p.u. 

The initial currents, voltages, flux linkages and torque 
angle are: 

i d0 = 0,286, i qQ = 0.640, v iQ = 0.384, v q0 = 0.97T, 
v f = 1.263, v 0 = 1.058, T f0 = 9.491, & 0 = 0.887. 
The time is scaled in computation hy a factor |3 = 7.308. 


numerical value 

of the 

system 

matrix 

A (eqn. 

(2.30)) 

is 

— 

0.0. 1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-0.683 -0.0346 

-0.0816 

0.0 

0.0 

0.882 

0.0 

1.32 

-0.0832 0*0 

-0.0254 

0.137 

0.0 

0.0 

0.0 

0.0 

0.130 0.0 

-0.107 

-0.137 

-1.37 

0.0 

0.0 

0.0 

0*0 0.0 

0.0 

0.0 

-0. 274 

0.0 

0.0 

0.0 

0.0 -0.0265 

0.0 

0.0 

0.0 

-0.0616 

-1.37 

0.0 

0.0 0.0 

0,0 

0.0 

0.0 

0.0 

-13.70 

0.0 

0.0 0.0531 

0.0 

0.0 

0.0 

0.123 

2.74 

-0.171 


(2.31) 


t 
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2.2.2 Power System Considered by Davison and Rau [21]: 


The power system considered in this section 
(Figure 2.4) differs from the system of Section 2.2.1. The 
synchronous generator is steam turbine driven and has one 
damper winding along each rotor circuit. The generator is 
connected to an infinite bus via a short line whose resis- 
tance is small and hence represented by an inductance only. 
Simplified representation for the governor turbine system 
is assumed. The control inputs considered are, u-^, the 
change in reference signal to the exciter and u^» the 
change in valve position. The state variables considered 
are, flux linkages of the field winding flux linkages 
of direct axis winding flux linkages of direct axis 
damper winding Y , flux linkages of quadrature axis 

winding ¥ , terminal voltage (measured) e. , field voltage 
q t 

E„, , mechitorcaie T , rotor speed w and the load angle 6. 
id ■ m 

Electrical transients in the electrical network, pl> p¥ 

d q 

and armature resistances are neglected. The performance 
of t£e machine is represented by eqns. (2.32) to (2.43). 



= x , i,,. + x , i, , - x, i, 
d ad fd ad kd d d 


w = x i, - x i 
M aq Hcq q q 

^ kd = X ad x fd + x kkd \d “ x ad x d 


(2.32) 

(2.33) 

(2.34) 

(2.35) 


( 2 . 36 ) 
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FIGURE 2.4 SYNCHRONOUS MACHINE INFINITE BUS SYSTEM 
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FIGURE 2.5 SIMPLIFIED EXCITER VOLTAGE REPRESENTATION 
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FIGURE 2.6 SIMPLIFIED VOLTAGE ERROR DETECTOR REPRESENTATION 
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mech. power3 input 
to the machine 


FIGURE 2.7 SIMPLIFIED POWER INPUT REPRESENTATION OF 

A THERMAL MACHINE 
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P kq " 

x, , i, —x i 
kkq kq aq q 

(2.37) 

> 

- “\a 

(2.38) 

kq = 

■’'kq ^kq 

(2.39) 

II 

<H 

CD 

¥ + wR i 

f d + id fd 

(2.40) 

II 

<H 

5H 

x ffd 1 fd + x ad *kd " x ad x d 

(2.41) 

• 

6 = 

to 

(2.42) 

o 

0* 


(2.43) 


On manupulat ing eqns. (2.32) to (2.43) (derivation omitted) 
the machine equations in the state variable form can be 
obtained . 

A sixth order system is required to represent the 
machine while each of the exciter voltage system (Figure 
2.5) » voltage error detector (Figure 2.S) and the governor 
turbine system (Figure 2.7) is represented by a first 
order system. Thus a ninth order model for the power 
system is obtained. The machine terminal conditions give in 
d-q- form, the following relation in linearized forms 


A I = [A i d A i^- 


and 


AI = 


-j v Sin6^ 


— — 

1- O 0 0 





/Jt * M 

d 

V v 


A¥ 

1 V Cos6^ 

0 1 ° ° 


4 

X x_ 


6 

e e 


' . t 


(2.44) 


The ninth order nonlinear equations representing 
machine, exciter voltage, voltage error detector and 
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governor turbine system are linearized around an operating 
condition and the linear model is obtained. For brevity 
the general form of the system matrix is not given here. 

The relevant data for the power system considered 
in this section are [2l]i 

Transmission Line 

Z e = 0.25 

Synchronous Machine (Typical values of a 500 MW set) 

\d = 1 * 89 ' x d = °- 27 ' x d = °- 175, T do = °*° 31 ’ 

T do = 4 °' T kd = 0 *°° 65 ’ ^d = °* 013 ’ T kdo = °* 392 ' 

x a<i = 1.89, x^ = 0.175, x fd = 0.128, R fd = 0.00124, 

x = 0.15, x q = 2.04, = 0.0261, x^ Q = 0.195, 

v Q = 1.0, H = 6.75, D = 2.36. 

Exciter and Voltage Regulator System 

K e = 155, = 0.01. 

Voltage Error Detector System 
x =0.01 

V 

Governor and Turbine System 

x = 0.3 sec, 

s 

The state vector x for the linear system is 

Z = [M w ,A, d . A, kd . A» q , Ae t , AB fd , SI,, A» . A 6] T 
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and controls are AE ~ and AT. _ • The linearized system 

rei mg 

obtained in standard state space form as 

k = A£ + Bu + Dv (2*45) 

where v is the disturbance vector. The numerical values 
of the matrices A, B and D are [2l]; 


B = D 


-3.07 

2.98 

0.0 

0.0 

0 

♦ 

H 

0.0 0.0 

-61.50 

0.52 

0.0 

0.0 

0.18 

0.0 - 0.09 

-15.50 

-32.20 

0.0 

0.0 

0.0 

0.0 0.0 

0.0 

0.0 - 

27.0 - 

•100.0 

0.0 

0.0 -0.08 

44.00 

0.0 - 

89.80 

-38.75 

0.0 

0.0 0.0 

0.0 

0.0 

0.0 

0.0- 

0 

0 

• 

0.0 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

- 3.3 0.0 

•223.00 

0.0 -47.80 

0.0 

0.0 

55.4 - 0.35 

0.0 

O 

* 

O 

0.0 

0.0 

0.0 

0.0 1.0 






(2 

0 0.0 

0.0 0. 

0 0.0 

25.0 

0.0 0 

1 T 

.0 0.0 

,0 0.0 

0.0 0. 

0 0.0 

0.0 

3.33 0 

.0 0.0 


— 0.6 

-32.1 

- 15.6 

-5.3 

0.0 

0.0 

0.0 

- 222.0 

o.j 

• 46) 


a = C 4 s ref - Ai ing ]- 


(2.47) 

(2.48) 


y = [ AE ref , AT . ng ]- 


2.3 DYNAMIC MODEL OF A TWO-AREA POWER SYSTEM [ 38 ] 

Linear model of a two-area power system (Figure 2.8) 
is used in this thesis to design the load frequency 
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controller. This model is the same as that considered by 
Fo sha and ELgerd [38] . The function of the load frequency 
controller is to keep the frequencies f-j_ and of the two 
areas and the tie line power flow P + . at the desired level 
for all load conditions and to restrict their variations 
within a close tolerance. Purt hermor e , the steady state 
value of the area control error in each area should be zero. 
This is achieved by controlling the inputs P c ^ and P^ of 
the generators. 

The linear model developed is based on the following 
assumptions: 

(i) There is no interaction between the load frequency 
control channel and the voltage control channel and 
the node voltages are treated as constant during the 
analysis . 

(ii) Strong electrical interconnection exists between 
generators of each area. Therefore, each area can 
be represented by an equivalent generator. 

(iii) The two areas are interconnected by a weak tie line. 

(iv) losses are neglected. 

(v) The system equations are linearized around an 
operating point. 

(vi) The turbine* governor and generator of each area are 
retire sehted by a first order system. 
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The state variables selected for the linearized 
model are 


x - [/AP tie , /A f-^, i£ lP AP gl ,A X gvl »/A f 2 , Af 2 t 


A P , A x ] 
g2 gv2 J 


T 


( 2 . 50 ) 


The change in tie line power in the linearized form 
can he represented as 


Ap = T ( 

tie 12^ 

/Af^ dt - /Af 2 dt ) 

(2.51) 

where 



m - 271 | 

T 12 x, 0 * 

Mil! cos(& 1 - a 2 | o 

(2.52) 


is a constant for the system defined as electrical stiff- 
ness of the line and is evaluated at the operating condi- 
tion. The differential equations representing the system 
(see Figure 2.3) performance are: 


pAP tie = T 12 ( dt - /a* 2 dt) (2.55) 

fS pi ^ = - [Hjl + 4P tle + 4E dl -A E gl ] (2.54) 

(2.55) 


pA p = - -i- AP , + ~~AZ , 

gl T +1 g 1 T +0 gvl 


pAX 


-tl 

1 


gvl 


■gvl 


■gi T t2 
A X 


■ f 4. ■■ —A P 

■s^ 1 VA 1 T 


'gvl 


cl 


p(/Af 1 ) = Af x 


(2.56) 

(2.57) 


where f is the nominal frequency. 
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For each, area 4 differential equations similar to 
eqns. ^2.54) to (2.57) are required to represent the dynamic 
performance of the system for small variation in load. Thus 
a ninth order linear model obtained for a two area powej? 
syst em is 

x = A x + B u (2.58) 

For brevity the general form of the system matrix is not 
reported here. . Also note that the origin of the state space 
of the above model is the steady state operating point. 

The relevant data for the power system considered 
in this section are [38]; 

Area 1 : 

Generator 

E ± = 5.0 , D 1 = 8.33 x 10“ 5 . 

Govern or -Turbine 

- °- 545 - Vi = °- 08> E i = 2 - 4 - 

Area 2 : 

Data for Area 2 are identical to those of Area 1. 
Tie-line; 

T 12 = 0,545 
Operating Condition 

< fi l - Vo - 30°- 



33 


The non-zero el$ments of the system matrix A and' 
control matrix B and. disturbance matrix x are: 

3 . 2 = 0*545, a^ g = -0. 545, a^ ^ = 1*0, a^ 2 ~ "“3*27 , 

a^ ^ = 0.05, a^ ^ = 6.0, a^ g — 3*27 , a^ ^ = —3*333, 


a 4,5 = 

3 . 333 , 

- “ 5 . 208 , 

p . — -12 • 5 9 

0 9 D 

a 6,7 = 

1 . 0 , 

a 7,2 = 

3 - 27 , 

a 7 , 6 ~ ” 3 * 27 , 

a rj rj » - 0,05 t 

a 7,8 = 

6 . 0 , 

a 8,8 = 

- 3.333 

* a 8,9 = 5 * 333 , 

a 9»7 * ~ 5 * 2 ' 

a 9,9 = 

“ 12 . 5 , 


b^ 1 — 12*5, 2 “ 12* 5 y 1 ~~~ ™6 . 0 and t y —6 • 0 * 

2.4 CONCLUSION 

In this chapter two models of power systems and 
one 2-area load frequency control model have been presented. 
These are used in subsequent chapters for illustrating the 
new modal control design techniques proposed in this thesis. 



CHAPTER 3 


OPTIMAL MODAL CONTROL OF POWER SYSTEMS 

3.1 INTRODUCTION 

In modern control theory there are mainly two 
approaches to the design of feedback controllers for 
linear, time-invariant dynamical systems to obtain satis- 
factory closed loop system behaviour: (1) procedures based 
on optimal control theory and (2) procedures based on pole 
assignment techniques. 

Optimal linear regulator theory for linear, time- 
invariant systems with quadratic performance indices has 
been widely used, to design feedback controllers for power 
systems [6-8,20-46]. In this approach the state and 
control weighting matrices in the performance index are 
selected iteratively and the resulting sequence of optimal 
control problems are solved to obtain satisfactory closed- 
loop behaviour. Satisfactory behaviour of: the closed loop 
system is checked either by studying the system response 
to certain test inputs [6-8,20,21,25], or by determining 
closed loop system poles [24]« The iterations involved in 
the optimal control theory approach are unattractive from 
the computational point of view. 

As in the classical root locus procedure, locations 
of the closed loop system poles give the designer a 'feel' 
for the behaviour of the system. Thus, direct design 
procedures which assign closed loop syst em poles to 
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desirable locations are to be preferred to the optimal 
control theory approach mentioned above. 

It is well known [66] that it is possible to achieve 
any arbitrary assigned closed loop eigenvalue locations 
(with the constraint that complex eigenvalues should occur 
in conjugate pairs), if the open loop system is controllable. 
In practice, however , the open loop eigenvalues which need 
to be shifted (that is, the dominant eigenvalues) are usually 
small in number. In such cases, it is not necessary that 
the system be completely controllable to obtain a feedback 
controller to effect the shift of the dominant eigenvalues: 
however, it is necessary that these eigenvalues be control- 
lable. Modal control theory [52-55] offers a useful tool 
for designing controllers for such systems. 

Modal analysis of dynamical systems is one of the 
best known and most extensively cultivated fields «>f 
classical physics. In the recent explosive development of 
control theory , the modal approach to the analysis and 
synthesis of multivariable control systems has led to many 
interesting and useful results [55]* The original adumbra- 
tion of modal control theory was given by Rosenbrock [67]. 
Following Rosenbrock 1 s work, many researchers developed 
various modal control procedures [52—55,68]. Some particular 
results relating to single input systems were obtained by 
ELlis and White [68]. The work of Simon and Mitt er [52] 
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constitutes a major contribution to modal control theory-. 
General closed form expressions for controller gains for 
single input systems were given by Porter et al [54-55] 
and Mayne et al [53]- 

In the power systems area, interest has been shp^m 
recently in direct pole assignment technique for controller 
design to obtain improved dynamic behaviour. Seraji [31] 
applied a pole assignment technique to a fifth order model 
of the synchronous generator system with a single input. 

Pai et al [49] considered a synchronous machine connected 
to an infinite bus with an improved representation of 
the exciter and the governor on the lines of Yu et al [6]. 
This resulted in an eighth order model of the power system 
with two control inputs. Applying modal control theory, they 
obtained feedback controllers and assigned dominant eigen- 
values to desired closed loop pole locations. Their pro- 
cedure involved utilizing only one of the available inputs 
to shift one or more eigenvalues at a time. This is a 
rather restrictive procedure and may result in high con- 
troller gains. 

In this chapter, considering a dyadic structure, 
for the feedback matrix an algorithm is presented for 
designing an optimal modal controller; that is, a 
controller in which the feedback gains have been mini- 
mized in a particular sense. Computational difficulties 
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and problems in designing controllers for large systems 
are also indicated. The algorithm developed is applied 
in designing control schemes for power systems, and the 
computational results are given. 

3.2 DEVELOPMENT 

Consider the linear, time-invariant, mult L- input , 

multi- out put system 
♦ 

x = A x + B u (3*1) 

Let the system eigenvalues X^» X^» **•> be distinct. 

It is desired to improve the system behaviour by shifting 

the 1 & 1 dominant eigenvalues, considered to be X X 

. . . , X^ f or convenience, by state feedback so that the 

closed loop system eigenvalues are (p^» Pg'**** P £ » 

X„ ,i»X ,o» *»*» X '}• It is assumed that the dominant 
x, +-L f XI 

modes are controllable [52]. The general problem of 
assigning the dominant poles to new locations without 
disturbing the locations of the rest of the poles and 
without constraining the structure of the feedback matrix 
P in any manner is a computationally difficult problem, 
especially in large systems. The problem is simplified 
by assuming that the feedback law is of the type [52], 

u=Ex=ag T x (3*2) 

where a and g are m and n-vectors respectively. The 
pole assignment problem for (3»1) can be seen to be 
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equivalent to the pole-assignment problem for the single 
input system 

x = A 2 . + £ 2 (3-3) 

where 

a = B a (3-4) 

and z = g T x (3*5) 


From (3*2) it can be seen that the elements 
a^, (i = l,2,...,m) of a denote the relative proportions 
in which the m control inputs u-^, u^» ...» u m are used, 
let a_^*s be so selected that the dominant modes of (3*3) 
are controllable [52]. It has been shown in Appendix A 

that for the desired pole assignment 

1 

g = I k. v. (3.6) 

- i=l 1 1 


where v.^, (i = 1,2,.,,, Jt ) are the reciprocal eigen- 
vectors of system matrix A associated with eigenvalues 
^1* ^2* * * • respectively^ the modal controller gains 
(i — 1,2,... ^ ) ar e 


k i ~ ^ 

0=1 J 


n (p, - 


K)M£ b a n (^ . 
1 1 0=1 3 


- x i)3 


(3.7) 


If we assume that in the closed loop system, complex 
eigenvalues appear in conjugate pairs, the vector g is 
then real [Appendix A] . 
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3.3 OPTIMAL MODAL CONTROL 

The feedback law (eqn.3.2) can be realized optimally 
so that the feedback gains, that is the elements of the 
modal controller matrix F are minimized. This may be done 
by choosing a^, (i = l,2,...,m) to minimize the performance 
ind ex 

m n 

j = I l CL 2 g? (3.8) 

X j=l i=l 3 1 

where the elements g^t (i = l,2,,..,n) are given by eqns. (3*6) 
and (3*7). Since for a given a and g, a g^=(r a)(^ g ? ), where 
r is any scalar, we may impose the constraints 

- 1 1 a j 1 + 1 > (0 = 1,2,. ..,m) (3.9) 

for computation convenience. Instead of choosing the 
objective function of eqn.(3»8) with constraints (eqn.3»9) 
we may choose the following simpler objective function for 
minimization without any constraints: 

m n 

J = j a* + 1 g\ (3.10) 

^ i=l 1 j=l 3 


Though the feedback gains are not directly mini- 
mized in the above objective function, it is clear from 
eqn.(3.10) that minimizing should yield small values 
for the feedback gains a g . , (j= 1,2, . . . ,m) , (i=l, 2,*. . , 
It is not necessary to impose any magnitude constraints 
such as in eqni(3.9). Any of the well known function 


n). 
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minimization techniques (see, for example [69])» may "be 
used to solve either of the minimization problems enunciated 
above. Computational results reported in this thesis 
have been obtained by using the gradient method described 
in Appendix B. The method of Grad and Brebner [70] based 
on the algorithm proposed by Wilkinson [71] is used for 
computing the eigenvalues and eigenvectors. 

3.4 ALGORITHM TO REALIZE THE OPTIMAL MODAL CONTROLLER 

i) Set iteration count n = 1. 

/ \ 

ii) Set objective function value ' = 0, where the 
objective function J is of the form (eqn. 3. 8 or 
eqn.3.10) and superscript n denotes the iteration 
count . 

iii) Determine the dominant eigenvalues and the corres- 
ponding reciprocal eigenvectors of the system matrix A. 

iv) Select new locations for the dominant eigenvalues. 

v) Generate £ from eqn. ( 3 . 4 ) by choosing a such that the 

T 

mode controllability condition, v. £ , 0, (3=1,2, ... ) 

j 

is satisfied. 

vi) Calculate k. , (j=l, 2, ...» & ) from eqn. (3.7) and 

J 

vector g from eqn. (3.6). 

( n+ 2 j 

vii) Calculate the new objective function value J v y 
from eqn. (3*8) or eqn. (3*10) as the case may be. 

viii) If | j( n+1 Lj( n ) | i s less than a specified small 
positive number* go to step (xi). 
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ix) Change a vector as a = a + c( A a ) where c is 

H cW 

step length and (Aa) is the negative of the 
gradient of the lagrangian function with respect 
to a [see Appendix B, eqns.(B.7) or (B.13)]. Check 
whether the mode controllability criterion for dominant 
modes is satisfied by the new a vector. If it is not 
satisfied, then change c appropriately to satisfy the 
criterion. 

x) Increment the iteration count n by 1. G-o to step (vi). 

xi) Calculate the modal controller matrix F from eqn. (3*2). 

xii) Stop. 

3.5 THE CONTROL STRUCTURE 

The general structure of the closed loop control for 
a multi- input, multi- out put (MIMO) system is shown in 
Figure 3.1. The stat e- variables are continuously monitored 
via sensors. The operation of the control system is based 
on signals which represent deviation of the states from the 
nominal settings. If the linear model is considered* then 
the variables to be considered are incremental quantities 
with respect to the operating point. Thus the models 
discussed in Chapter 2, Sections 2.2 and 2.3 can be inter- 
preted in terms of the block diagram (Figure 3»1) appro- 
priately. To illustrate application of the optimal modal 
controller algorithm -O? Secs.3 . 3 4 , two power systems 
problems are considered in this chapter: 



Disturbances 
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(a) Controller design to improve dynamic response of 
a power system, and 

(b) load frequency control problem for a two-area 
power system. 

3-6 DESIGN OF MODAL CONTROLLERS FOR POWER SYSTEMS 

In this section, we shall apply the theory of 
Sections 3*3 and 3*4 to improve the dynamic performance 
of a power system. Two systems whose linear mod-els were 
presented in Chapter 2 are considered. 

(1) Synchronous machine connected to an infinite bus 
and driven by hydraulic turbine. 

(2) Multi-area load frequency system. 

3.6.1 Synchronous Machine-Infinite Bus system: 

The power system considered consists of a synchro- 
nous generator connected to an infinite bus through along 
transmission line. The generator is driven by a hydraulic 
turbine. The system has an exciter voltage regulator and 
a turbine speed governor. The generator dynamics are 
represented by a third order system, while the exciter 
voltage regulator and primemover-speed governor dynamics 
are represented by second order and third order system 
respectively. This power system was considered by Yu et al 
[ 6 r ] for designing an opt ima-L regulator while Pax, et al [49] 
.designed a modal controller for this system. <. The mathema- 
tical. model, of the system has been, discussed in Section 
2.2.1. The state space model is given by 
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z = A x + B u 

where A and B are given hy eqns. (2.28) and (2-30) respec- 
tively. The operating points correspond to 0.735 p.u. real 
power, 0.034 p.u. reactive power delivered by the machine 
and 1.05 p.u. terminal voltage. The eigenvalues of the 
uncontrolled system matrix A are 

-0.0114 ± 3 0.7986; -0.0572; -0.0772+3 0.1146; 

-0.1952; -0.2740; - 13.70. 

From these eigenvalues it can be inferred that 
though the uncontrolled system is stable, the response 
will be unsatisfactory due to the locations of the first 
five eigenvalues very near to the 3 w axis. These eigen- 
values are termed as 'Dominant eigenvalues'. The uncon- 
trolled system requires considerable time to settle down 
to a steady state value when a disturbance takes place 
(see figure 3*2). The disturbance considered is the 
release of a load perturbation A P which gives the initial 
conditions [49] as 

z(o) =[ 0.0258, 0.0, 0.0092, 0.0171, 0.0, 0.0, 0.0, 0.0 ] T 

The performance of the system is proposed to be improved 
by designing a modal control to shift the dominant eigen*^ 
values to the following locations* 

-0.40±3 0.915; -0.60; -0. 96+3 0.720. 



.Torque angle in radians A 6 


0.02 1 




a-. Without controller 

h _ With optimal modal controller 
case I 

c - With optimal modal controller 
case II 


FIGURE 3,2 SYSTEM RESPONSE FOR A <5 
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The locations of these eigenvalues has been justified in 
[49] on the basis of sector criterion. 


The optimal modal controller algorithm developed 
in Section 3*4 is applied to design a modal controller 
for the above syst em. Lagrange multiplier formulation and 
gradient technique are used for minimizing the cost func- 
tion. The resulting modal controller matrices are as 
follows; 

Case I: 



-1.384 -2.612 4.132 -6.691 58.44 -209-5 16. 31 -23.46 

5-221 9.853 -15.585 25.233 -220.45 790.3 -61.53 88.50 


when the cost function is optimized (eqn.3.8). 


Case II ; 

f 


F = 


-1.396 -2.619 4-138 -6.696 58.49 -210.7 16. 41 -23*53 

5.237 9.825 -15.521 25.120 -219.42 790.3 -61.58 88.27 


when the cost function is optimized (eqn. 3. 10). 


as 


The norm of the above feedback matrices defined 
m n 


J E f?. are 734019.63 and 734017.94 respectively. 
i=l j=l 13 


The transient response characteristics of the power 
system with optimal controllers obtained are shown in '■ 
Figure 3*2* The disturbance considered is the same as for 
the uncontrolled system. Standard 4th order Runge- Eutta 
method was used for solving the differential equations of 
the system. It is noted that the gain matrix F for both 
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the cases are almost identical. However, computationally 
minimization of is preferred since lias a simpler form. 

3.6.2 load Frequency Control of a Two -Area LFC System: 

The two area LFC system considered is the same as 
the one discussed "by Elgerd and Fosha [38]* The two areas 
which are considered identical are connected via a tie-line. 
Each area is represented by an equivalent generator. A 
load frequency controller is to be designed to keep the 
frequencies in the two areas and the tie-line power flow 
at the desired level for all loads and to minimize their 
deviation during a disturbance. This is done by control- 
ling the inputs to the machines in the respective area. 

The development of the 9th order state space model for the 
system has been presented in Section 2.3 where the numerical 
data are also given. The state space model is of the form 
x = A x + B u . 

The eigenvalues of the system matrix A are 

0.0; 0.0; -0.4972+03.522; -1.296+02.512; -1.624; 

-13.26; -13.29. f 

The "uncontrolled system is stable but requires 
«^i?trtime to stabilize when a step load disturbance of 
0.01 p.u. takes place in area one [37, 3-8*42] # The 

performance of the system is to be improved by designing 
a modal controller to shift the dominant eigenvalues to 
# the left, away from the j a -axis. 
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The set of dominant eigenvalues contains a pair 
of repetitive eigenvalues located at the origin. The 
algorithm developed in Section 3.5 cannot he applied to 
shift repetitive eigenvalues. This difficulty can he 
overcome either hy having a reduced state space model 
[ 42 ] or hy providing a suitable feedback which yields a 
set of dominant, distinct eigenvalues. This feedback 
should he such that the locations of the nondominant 
eigenvalues are not disturbed appreciably. The algorithm 
of Section 3.5 can then be applied to design a modal 
controller for this modified system. The resultant 
controller matrix is obtained by summing the controller 
feedback matrix for the modified system and the feedback 
matrix used for converting the repetitive eigenvalues 
into distinct eigenvalues. The latter approach has been 
used here in designing the load frequency controller. 

Considering the first order eigenvalue sensitivity 
[56,72] and the structure of control matrix, it is found 
that the following feedback matrix 

0. -0.96 0 0 0 0 0 0 0~j 

0.01 0 0 0 0 0 0 0 Oj 

shifts the two repetitive eigenvalues of zero magnitude 
to two distinct locations. The eigenvalues of the 
modified system after this feedback are 

-0.01017} ~0.1066iJ3.5l2| *-0.7537* -0. 76 2l±;j 2.531} 
-2.832} -13.15 I -13.28. 




49 


The system matrix of the modified system is 
I = A + B ? 1 

Now an optimal modal controller feedback matrix is 

designed so that the eigenvalues of the resultant closed 
loop system are 

-0.7537, -2.0, -2.18+34.201, -2.7±33.629, 

-2.832, -13.15, -13.28. 

The locations of the two pairs of complex poles are the 
same as that obtained by Ramamoorty et al [ 42 ] and the 
location -2 is arbitrarily selected so that it lies well 
within the set of nondominant eigenvalues of matrix A. ^he 
other locations are the same as those obtained for mat rix A. 

The algorithm of Section 3*4 is then applied to 
design an optimal modal controller to obtain the above 
closed loop spectrum. The optimal modal controller 
matrix is 

Case I i Minimization of J-^. 

- 3.0 - 2.22 0.26 1.73 0.52 3.18 1.48 1.69 0.38 

F 2 = 

* 10.53 7.75 -0.89 -6.05 -1.81 -11.10 -5.17 -5.91 -1.32 

and the resultant modal controller matrix is 


F 


-3.00 -3.18 0.26 1.73 0.52 3.18 1.48 1.69 0.38 


£10.54 7-75 -0.89 -6.05 -1.81 -11.10 -5.17 “5*91 -1.32 

m n 

The norm of feedback matrix F defined as I If?, is 435* 

i=l 3=1 13 
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Case II : Minimization of 

-0.62 -0.55 0.09 0 44 0.15 0.7^ 

P = 

2 10.46 9*05 -1.55 -7.26 -2.12 -12.78 


0.29 0.27 0.56 

-5.01 -4.61 -0.93 


and the resultant feedback matrix is 


F = 


-0.62 -1.51 0.09 0.44 0.13 0.77 0.29 0.27 0.56 
10.47 9 .05 -1.55 - 7.26 -2.12 -12.78 -5-01 -4.61 -0.93 


The norm of feedback matrix is 463 . 

The transient responses for the controlled system 
with above modal controllers are shOT-m in Figures 3»3 and 
3»4 for a step load disturbance of 0.01 p.u. in area one. 
The transient responses compare favourably with the respon- 
ses given in reference [ 42 ] for the optimal controller. It 
is also noted that the two feedback matrices obtained above 
have nearly the same matrix norm. 

3.7 CONCLUSION 

A simple method for designing regulators based on 
modal control theory has been suggested to obtain improved 
performance of linear time-invariant systems. An efficient 
and straightforward algorithm to design optimal modal 
controller has been developed. The algorithm is applied 
to design suitable controllers for (a) improving dynamic 
response of a power system and (b) load frequency control 
of a two area system. The algorithm has also been applied 
to shift repetitive eigenvalues in the dominant set in the 
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case of load frequency control problem. In this case, 
the repetitive eigenvalues are converted to distinct ones 
by a suitable feedback in order to enable application of 
the modal control algorithm. 

Computational results show that in the problem 
considered, large feedback gains are required to shift 
all the dominant eigenvalues at one time to the desired 
closed loop locations, large feedback gains are undesir- 
able due to equipment limitations and noise amplification 
in the feedback path. It is, therefore, desirable to 
explore the possibilities of designing modal controllers 
with smaller gains. The problem of having reduced gains 
assumes more importance in designing controllers for large 
systems, where the number of dominant eigenvalues may be 
large. The dyadic structure for the controller matrix 
assumed in this chapter is rather too restrictive 
especially in high order systems with large number of 
inputs. In the next chapter we will discuss some simpli- 
fied techniques for large dynamical systems. 



CHAPTER 4 


OPTIMAL MODAL CONTROLLERS FOR LARCE SISTEMS 
4.1 INTRODUCTION 

The power utilities are generally interconnected to 
form large grids for economic and reliable operation. These 
interconnections result in large dimensions of the state and 
control vectors of the system models. From the computational 
point of view, the conventional feedback design techniques 
may not be suitable for these high order systems. Davison 
and Chadha [73] have established the advantages of the modal 
control approach for large systems. 

The problem of developing suitable techniques for ..... 
the design of controllers for large systems is an important 
and difficult one. Schemes based on modal reduction technique 
have been suggested for this problem [42, 74-76]. However, 
the closed loop system performance seems to be excessively 
sensitive to the reduced order model chosen and it may even 
be difficult to ensure satisfactory performance of the closed 
loop system [73» 74, 77]. It is- desirable that, in addition 
to the technique developed being computationally attractive, 
it should also provide the designer with an insight into the 
problem which can enable him to utilize his -engineering 
judgement to develop designs which are suitable from the 
practical view point. Modal control approach has this 
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feature [75] • Even though the approach of Davison and 
Chadha [73] is suitable for high order systems, the structure 
assumed for the controller is highly restrictive. This may 
result in rather large feedback gains. The same limitation 
exists in the approach given in Section 3.4* In this 
chapter, methods will be developed to overcome these limitation! 

The simplicity of the approach of modal controller 
design developed in chapter 3 is essentially due to the dyadic 
structure assumed for the feedback matrix. This restriction 
on the structure of the feedback matrix also permits a simple 
procedure for optimizing feedback gains. As was mentioned 
in chapter 3, the general problem of determining optimal modal 
controllers in which no constraint is imposed on the structure 
of the feedback matrix is a computationally difficult one, 
especially for large systems. In the technique developed in 

this chapter, the modal controller feedback matrix is obtained ; 

[ 

as a sum of dyadic matrices. Thus the structure of the 
feedback matrix is more general than those assumed in chapter 3i 
This structure retains the simplicity and computational ! 

advantages of cthe dyadic structure in feedback optimization I 

I 

and yields lower feedback gains. It permits the designer to j 
use his engineering judgement and intuition in developing ! 

suitable designs. In many situations, all the inputs need | 

I 

not have to be used for feedback control: some : of them may 
be discarded without affecting the feedback gains appreciably. 
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The procedures of this chapter help the designer to identify 
the inputs to he used for feedback in such situations. 

In essence, in the techniques reported in this 
chapter, the dominant eigenvalues are divided into groups 
and the algorithms of chapter 3 are applied sequentially to 
the various groups to obtain the required locations of closed* 
loop poles. Thus for the pole allocation problem for any of 
the groups, the < dyadic structure of the feedback matrix is 
retained* In many situations, dividing dominant eigenvalues 
into a s mal 1 number of groups will yield structures significantly 
better than those obtained without grouping. It can thus be 
seen that algorithms being proposed are slightly saore complex 
than those of chapter 3* 

4.2 DEVELOPMENT OP METHOD BASED ON GROUPING 
4.2,1 Eigenvalue Assignment in Groups 

Consider the system defined by eqn, (3*1)* bet the 
•A' do minan t eigenvalues of the uncontrolled system be 
divided into *p’ groups as follows: 1 

Group 1 : * ^ 2 ^* **•» * Al"*^ ^ 

Group .2 : ^ 

• * 

• * 

Group p : » ^2^^ **•» *j,p^ ^ 
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where $, is the number of eigenvalues in q^ group. 

Then, clearly H=Z + & + 9 . . i 

J. d p 

It is assumed that each group has distinct eigenvalues. Let 
the corresponding desired locations for the closed loop poles be 


Group 1 : [ , P^ 1 ), ..., p £1 ^ 

Group 2 : [p-^ 2 ) , P 2 ^, 


] 

] 


Group p : [p 1 ^ p \p 2 ^ p) , ...,p Ap (p) ] 

In any of the groups above, complex eigenvalues appear in 
conjugate pairs. 


Wsing eqn. (3.2 ) e a dyadic matrix 

p(D =a (l) g (l) T 


,' 1 * 

(4.1) 


to shift the eigenvalues of Group 1 to their desired locations 
can be obtained. The resulting closed loop system matrix 
after shifting the eigenvalues in Group 1 is. 


A 


( 1 ) 


= A + B 

= A + B g^ 


T 


(4.2) 


Similarly a feedback matrix 

I<2> = a< 2 > g( 2 > T 

to shift the eigenvalues of Group 2 to their desired locations 
can be computed. However, in the computation, the reciprocal 
eigenvectors of the new system matrix A^^ given In eqn. (4.2), 
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corresponding to the eigenvalues of Group 2, are to he 
used. Thus in general, to shift the eigenvalues in the 
Group q after the assignment of the eigenvalues of the first 
(q- 1 ) groups, the required feedback matrix is 

p (q) = g(q) T f ( q = 1 , 2 ,..., p) (4.3) 

where 

j k.^ v.^ , ( 4 . 4 ) 

i=l 1 1 

k i (4) = f jii < Pj (q) : x i <4) )J / C2i t,l)T B a (9) 

1 < x d (4) - x i W)] > 

(i=l, ..., ft q ) (4.5) 

(i=l, 2 ) .... ft^) denote the reciprocal eigenvectors 
of j^G”- 1 -) associated with eigenvalues (i=l, 2 > . . . , ft ^) 

and the matrix A^^ is 


A<»> 

+ 

«! 

1! 

q. 

i 


(4.6) 



±=i 

mom 


The elements 

(i=l,2. 

« » • i m) of a. 

denote the relative 


proportion in ’which the m control inputs are used in shifting 
the eigenvalues in the q group. These elements are to be 
so selected that the mode controllability condition 

v.^ T B 5 ^ 0 , (i=l, 2 , ..., & ) is satisfied 

—1 — q 
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The resultant feedback matrix F is the sum of the 
p dyadic feedback matrices. Thus 

F (i) 


P 

E 

i=l 


P (i) (i) T 
E a K ' g K J 

i=l 


(4.7) 


The rank of F is less than or equal to m. 

The feedback law (4.7) can be realized optimally so 
that the feedback gains , that is elements of the modal 
controller matrix F are minimized. This may be done by 


(4.8) 


choosing 

a. (4) , 
1 ' 

(i=l,2j..,m) 


■ n 

m , P f, 

J 

E 

i=l 

E ( X. 

3=1 *- L 3 


The Lagrange multiplier technique may be used to solve the 
optimization problem as in Chapter 5. 


4.2.2 Grouping Criterion 

A criterion- for grouping of dominant eigenvalues may 
be obtained from an examination of the equation which yields 
feedback gains for the case where the eigenvalues are not 
grouped [Chapter 3* eqns.(3.2), (3.6) and (3.7)]. ^he equations 
are reproduced below for convenience: 


u = F x = a g x 


g = X K Z ± 
i=l 1 1 

a 


r - V] ' tJi T 


(5.2) 

(5.6) 

l 

& (X.-X )] (3.7) 

3-1 D - 

3f 1 (i=l,2, A ) 
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For given iiOGations of eigenvalues, to obtain low feedback 
T 

gains k^, v^ B a should be large for i=l,2,..., £ as seen 
from eqn. (3.7). In other words, the vector a should be 
so chosen that the inner products of the vector (Ba) with 
Xj_> (i=l,2, . . . , & ) should be large. There is obviously a 
limitation to the gain reduction that can be achieved by 
arranging choice of a in this manner. 

Suppose now we group together those reciprocal 
eigenvectors (assumed to be normalized) such that 

m 

product v. v . is close to unity for -each .of 
the p groups. Then if the eigenvalues corresponding to any 
one of these p groups are moved by the algorithm of Chapter 3# 
’fehe corresponding terms v^ B a will be uniformly large and 
the feedback gains correspondingly small. It is clear that, 

m 

for any group, the vectors B v^ will be 'nearly* collinear. 

From the arguments given above, it is clear that the 
mode controllability matrix 

M = B , (4.9) 

where the columns of the matrix are the reciprocal 
eigenvectors of the system matrix A corresponding to the 
dominant eigenvalues, provides a suitable basis for grouping 
the eigenvalues. The eigenvalues corresponding to the ___ 
’nearly* linearly dependent rows of M are grouped together. 
For considering the linear dependency of the rows of matrix 
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M, the elements of the modal matrix M are replaced by their 
absolute values. The method is explained in terms of the 
numerical example in Section 4.4. 

The feedback gains for the modal controller matrix 
will be generally less as the number of groups is increased. 
But the computational effort increases also with the 
increase in the number of groups. Therefore a compromise 
is to be made between the computational effort and the 
feedback gains while deciding upon the number of groups. 

4.3 ALGORITHM TO REALIZE THE OPTIMAL MODAL CONTROLLER 
FOR LARGE SYSTEMS 

We outline below the algorithm to move eigenvalues 
based on grouping technique. 

i) Determine the dominant eigenvalues and the 

corresponding reciprocal eigenvectors of the 
system matrix A. 

ii) Obtain the mode controllability matrix from 

eqn. (4.9). Divide the eigenvalues to be shifted 
into p groups using the criterion of Section 4.2.2* 
Arbitrarily order the groups. 

iii) Let n denote the iteration count in the optimization 
process. Set n equal to 1. 

iv) Set objective function = 0. 

v) Choose initial value of aV^ such that the modes 
of the first group are controllable. 
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vi) 

vii ) 

viii) 

±x) 

x) 

xi) 

xii) 

-Mi) 

xiv) 

xv ) 
xvi) 


Set q equal to 1. * g' denotes the group whose 

eigenvalues are being assigned. 

Compute from eqn. (4.3) through eqn.(4.4) and 

(4.5). 

Obtain the closed loop system matrix A^^ from 
eqn. (4.6 ) . 

If q equal to p, then go to stepfxiii),* 

If n=l, choose the initial values of such 

that the modes of (q+l)" 1 '' 11 group are controllable. 

Go to step (xii) • 

Check whether the mode controllability criterion for 
dominant ;modes in (q+l) group is satisfied. If 
the criterion is not satisfied then go to step(xvii) « 
Increment the iteration count q by 1 and go to 


setp vii. 

Calculate the new objective function j( n+1 ) from 
eqn. (4.8). If 1 J vn+1 ^- I is less than a 

specified small positive number, then go to step(xviii) . 
Change vector as 2^^/^ = 2^ + c ( A 2^)» 

( i=l ,2,...,p), where c is the step length and 
( Aa^) is the negative gradient of the Lagrangian 
function with respect to 
Increment the iteration count n by 1. 

Check, whether the eigenvalues of first group are 
controllable. If controllable then go to step(vi). 
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xvii) 


xviii) 
xix ) 


Change the step length c and compute = a , 

/ . \ \J-LU. 

c (ia 1 ), ( i=l , 2 ; . . . , p ) . The value of c should he 
such that the modes of first group are controllable. 
Go to step(vi) . 

Obtain the modal controller matrix P from eqn. (4.7). 
Stop. 


+ 


4.4 DESIGN OF OPTIMAL MODAL CONTROLLERS FOR POWER 
SYSTEMS BY EIGENVALUE GROUPING 

Two examples of the design of optimal, modal 

controllers by eigenvalue grouping are given in this section. 

One of them relates to the problem of improving the dyn am ic 

response of a power system, while the other pertains to the 

two area load frequency problem. 

4.4.1 Example 1: Power System 

The power system considered is the same as in Section 
3.6.1. The desired closed loop pole locations are the same 
as those considered in Section 3.6,1. 

The mode controllability matrix defined by eqn. (4.9) is 


-1.7xl0” 2 - o5.2xl0“ 3 

1.32X10" 1 

+ 3l.67xl0“ 2 

-1.7xl0” 2 + 3‘5.2xlO“ 3 

1.32xl0~ 1 

- jl.67xl0” 2 

-l.lxlO” 2 

0.96X10 -1 


0.99+30.0. 

-5.0x10” 2 

+ j4.05xl0“ 3 

0.99-30,0 

-5 . 0x10 ” 2 

- j4.05xl0” 3 
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It is observed that 


M ul l M ail * Ki 


and 


K 2 I 

K 2 1 

l M 4l| 


M 


22 


M* 0 t 

32 i 


M 




M 


51 


.th 


where m. - is the element of matrix M corresponding to i row 
and column. Therefore rows 1,2 and 3 may be considered 
nearly linearly d ep end ent.If sing the grouping criteria of 
Section 4.2.2, the dominant eigenvalues are divided into the 
two following groups. 


Group I : [-0.0115 ± jO.915, -0.60] 

Group II : [-0.96 + J0.720] 

The algorithm detailed in Section 4.3 is used and the 
cost function defined in eqn. (4.8) is optimized, She optimal 
modal control feedback matrix F is computed to be 


-0.474 -0.733 1.910 1,312 -1.853 -1.996 -0.053 -1.289 

F = 

0.294 -5.250 2.094 1.296 -1.794 5.328 -1.232 -3.765_ 

Let the norm of the matrix F = £ s’ f ,• ^ (4.10) 

i=l 3=1 


where f . . is the element of matrix F corresponding to i 
row and column. 


.th 


The norm of above matrix F is 96.22. Thus it can be 
seen that considerable reduction in feedback gains is obtained 
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by assigning the eigenvalues in groups (see col. 1,2 and 3 
of Table 4.1). The response of the power system with this 
modal controller matrix is better than the modal controller 
feedback matrix obtained by shifting all the eigenvalues in 
one group (see Fig. 4.1 and Fig. 3.2). 


4.4.2 Example 2 : l.F.G Problem 

A two area power system is considered for designing 
modal controller for load frequency control of the system. 
This problem was also considered in Section 3.6.2. 

The dominant eigenvalues are divided into two groups 
as mentioned below. 


Group I : [-0.102, -0.1066 + 33.51 2 ] 

Group II : [-0.762 +J2.5310] 


The modal controller feedback matrix F, obtained by- 
utilizing the algorithm of Section 4.3 and considering the 
cost function defined in eqn. (4.8) is 


F = 


-0.468 -2.477 -1.960 -2.024 -0.406 -0.749 -0.089 
5.204 1.245 1.154 0.861 0.132 -1.226 -1.054 


0.304 0.112 


-1.604 -0.396 


The norm of the feedback matrix F as defined in 
eqn. (4.10) is 29.57* 



with Optimal modal controller Sec. 4 .4.1 
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The feedback gain matrix when eigenvalues are 
shifted in groups has a smaller norm than that obtained 
for the case when eigenvalues are shifted in one group only. 
Response of the closed loop system with the initial condition 
of Section 3.6.2 is given in Rig. 4. 2 and ^ig.4.3. 

4.5 SUB-OPTIMAL MODAL CONTROLLER 

The algorithm detailed in the previous section 
yields small feedback gains for large values of p where p 
denotes the number of groups. However, more computational 
effort will be required in solving the optimization problem 
because in the iteration process for shifting the eigenvalues 
in any group, the reciprocal eigenvectors associated with 
the eigenvalues in that group are to be computed. Therefore 
for large systems, where the dominant eigenvalues may have 
to be divided into large number of groups, the above procedure 
may not be computationally attractive. For such systems a 
procedure of designing sub-optimal modal controllers which 
needs smaller computational effort is presented in this 
section. 

The dominant eigenvalues are divided into p groups. 
Instead of optimizing the resultant feedback matrix F after 
shifting all the eigenvalues, the optimization procedure is 
carried out in stages i.e. the feedback matrix required to 
shift the eigenvalues in each group is optimized. The resultant 
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FIGURE 4.2 SYSTEM RESPONSE FOR Af«j 
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feedback matrix is determined by summing the p feedback 
matrices thus obtained, that is, 

F = z P, (4.11) 

i=l 1 

where is the optimum feedback matrix to shift the eigenvalues 
of i group to the desired locations with eigenvalues of the 
first £i-l) groups already assigned. In this approach the 
reciprocal eigenvectors are computed only for a total of p 
times. Either of the performance indices defined in eqns.(3.8) 
and (3.10) may be selected for the optimization problem. 

4.6 DESIGN OF SUBOPTIMAL MODAL CONTROLLERS FOR 
A POWER SYSTEM 

Two examples of the design of subopt imal modal 
controllers are considered in this section to demonstrate the 
simplicity of the method explained in the above section. 

4 . 56.1 Example 1 

The system of Section 4.4.1 is considered for 
designing suboptimal modal controller. The distribution of 
eigenvalues in groups and their desired closed loop locations 
are identical to those in Section 4.4.1. The feedback matrix 
F defined in eqn.(y»^ll) is computed as detailed above. 

COMPUTATIONAL RESULTS 
Case I: Optimization of (eqn. (3.8)) 
a (1) = [-0.4501 4,118]* 

a (2) = [ 4.085 0,5714]* 
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F 


1.297 0.998 1.924 1.226 -1.707 -10.920 0.224 -4.2?4 
1.053 -4.895 0.353 0.071 -0.068 6.689 -1.289 -3.253 


The norm of the matrix F defined by eqn. (4.10) is 233*23 


II: When the cost function defined in eqn. (3.10) 
is optimized: 

a (1) = [-0.3542 3.176] T 

a (2) = [ 3.384 0.4562] 1 

1.314 1.014 1.927 1.227 -1.780 -11.015 0.226 -4.266 

F = 

1.048 -4.895 0.344 0.065 -0.059 6 ? 928 -1.289 -3.234 


The norm of the matrix F defined by eqn. (4.10) is 

236.07. 

The response of the closed loop control system with 
the sub optimal controllers design (see Fig. 4.1) are comparable 
with the system response with optimal modal controller (Section 
4.4.1) for identical disturbance and initial conditions. The 


norm of feedback matrices of the sub optimal designs is more 
than that of the optimal design. 


4.6.2 Example 2: 

The two area power system of Section 4.4.2 is 
considered for designing the suboptimal load frequency 
controller. The distribution of eigenvalues in groups and 
their desired closed loop locations are the same as in 
Section 4.4.2. 
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The cost function ^ defined in eqn.(3.10) is 
minimized while computing the feedback matrix for each 
group. The resultant suboptimal modal controller feedback 
matrix is 


P 


3.5079 -1.4410 -3.3881 -3.9321 -0.826 -3.625' 1.4417 
3.8534 0.3214 -1.1761 -1.5190 -0.334 -2.9036 -0.1427 

1.7692 0.3481 

0.1047 0.29. . 

The norm of tthe matrix P defined by eqn. (4.10) is 88.83. 


The response of the closed loop control system with 
suboptimal controller (see Pig. 4. 2 and 4.3) is comparable with 
the system response with optimal modal controller (Section 4.4*2] 
The norm of the suboptimal feedback matrix is higher than that ; 
of optimal feedback matrix. 

4.7 NEAR OPTIMAL MODAL CONTROL WITH DOMINANT CONTROL INPUTS 


Inspite of the improved technique suggested in | 

Section 4.6 computational effort may still be considerable. j 

This is due to the large dimension of the a vector associated 
with each group of dominant eigenvalues with respect to 

} 

t : 

which the objective function has to be minimized f However, > 
the computational effort may be reduced by associating 

with each group of dominant eigenvalues a corresponding set I 

' ... • ' [ 

consisting of a small number of inputs to be used for j 

| 

shifting the eigenvalues in the groups optimally. When 
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these inputs are properly selected, the objective function 

is marginally effected. Such inputs are termed as 'Dominant 

Control inputs’. These input sets may be formed by studying 

the 'degree of coupling' between the modes in a given group 

and the inputs. The degree of coupling between the i mode 
th 

and j input may be considered to be given by 


Si 


T 




where v^ is the normalized reciprocal eigenvectors of A 
associated with X^ and b.. is the j column of matrix B, 

Let M be the mode controllability matrix associated with 
the q group, that is, 

M = V T B 

q. a 

where the columns of are the normalized eigenvectors 
associated with the group. Then the elements of M^ indicate 
the degree of coupling of the eigenvalues in the group with 
the inputs. The inputs corresponding to the columns of M^ 
having- relatively large magnitudes form the set of dominant 
inputs. 


4.8 NEAR OPTIMAL MODAL CONTROL OP POWER SYSTEMS 

The examples of the design of near optimal modal 
controllers are solved in this section to demonstrate the 
suitability of tfce method suggested in Section 4.7. 
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4.8.1 Example 1 : 

The example of Section 4.4.1 is considered here with 
the same eigenvalues grouping. 


The mode controllability matrices for the two groups are: 



-1.7xl0~ 2 -j5.2xl0“ 3 
-1.7xl0“ 2 +j5.2xlO“^ 
-l.lxlCT 2 


1.32X1CT 1 + ji.67xl0" 2 
1.32X10" 1 - jl.67xl0“ 2 
0.96xl0"' 1 



0.99+j0.0 

0. 99-jO. 0 


-5.0xl0~ 2 +j4.5xl0“ 5 J 
-5 . OxlO" 2 - j4 . 5xlO“ 3 J 


On the basis of the degree of coupling (Section 4.7) between the 
modes in a group and the inputs it is found that the dominant 
control for Group I is and the dominant control for Group II j 
is u-^. Since for shifting eigenvalues of each group only one 
dominant control, input exists, there is a unique resultant 
feedback matrix for the pole assignment problem. The near-optimal' 

i 

modal controller feedback matrix is found to be 


P 


0.481 0,184 ’1.818 1.204 -1.766 -6.566 -0.093 -2.774 
0.881 -5.028 0.086 -0.100 0.183 8.380 -1.320 -2.678 


The norm of the feedback matrix P as defined by 
eqn, (4.10) is 164.2. It is better than that of suboptimal 
design and is near to the value of the optimal feedback 
matrix (Section 4.4.1). Response of the closed loop system 
with near optimal modal controller for the disturbance 
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considered in Section 4.4.1 is plotted in figure 4.1. 

4.8.2 Example 2: 

In this example the dynamic stability of steam turbine 
driven synchronous generator is improved by designing a near 
optimal modal controller. 

The power system considered consists of a steam 
turbine driven synchronous generator connected to an infinite 
bus through a radial transmission line. The synchronous 
generator has damper winding along the direct as well as 
quadrature axes. Both the exciter and governor inputs are 
available as control inputs. This system was considered by 
Davison and Rau [21] for designing an optimal controller. 

The system and control matrices for the linearised model and 
the numerical values of parameters are given in Chapter 2 
Section 2.2.2. The eigenvalues of the system matrix are 
-0.101 + jlO.47, -3.35, -13.55, -29.62, -33.867 ± jlO.239, 
-106,855 + jl5.889 and the desired dominant eigenvalues "... 
of the closed loop system are specified to be -9.07 + 0*17.3, 
—10.8, -17.3. These locations are the same as obtained by 
Davison and Rau [21]. The dominant eigenvalues of the 
uncontrolled system are divided into following two groups. 

Group I : [-0.1010 + jl0.47, -3.33] and 

Group II : [-13.55] 
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It is found that the dominant inputs for Croup I and 
II are u^ and u^ respectively. The resultant feedback 
matrix of the modal controller is found to be 


[-14.2 6.4 -2.04 -42.24 6.5 -0.149 1.6 0.44 

F = I 

j _ 8.7 7.5 0.36 2.41 -2.4 0.081 -7.6 -2.50 

The feedback matrix obtained by Davison and Rau [ 21] 


4.66 

-4.74 

is 


-1.45 0.14 -0.26 -1.94 2.62 -0.14 0.52 0.005 0.41 


F = 


[11.40 7.62 0.60 4.50 -1.98 0.07 -8.64 -2.770 -9.44J 

It appears from the above that Davison and Rau’s 
result is superior. However it must be noted that their 
feedback disturbs the position of the nondominant eigenvalues 
to such locations as to give better feedback gains. Responses 
of the closed loop system with near optimal modal controller 
are plotted in Figures 4*4 and 4.5. 


4.9 AIT ALTERNATE METHOD OF MODAL CONTROLLER DESI&N 
FOR LARG-E SYSTEMS 

Reciprocal eigenvectors are to be computed (p-l) 
times in each iteration of the optimization process described 
in Section 4»5 with the cost function J defined in eqn. (4.8), 
where p is the number of groups in which the dominant 
eigenvalues are grouped. This may require more computer time 
for each iteration, especially in the case of large systems. 
Therefore to decrease the computational time, either the 
number of groups should be small or suboptimal modal 
controller design procedure (Section 4.6) is to be used. 



Torque angle in radians A 6 
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Both the approaches result in increased feedback gains. To 
overcome the above difficulties, an alternate procedure of 
designing optimal modal controller which retains the 
computational simplicity of subopt imal design procedure 
(Section 4.6) and permits the division of the eigenvalues 
into large number of groups to yield better feedback gains 
is presented in this section. In this procedure reciprocal 
eigenvectors are to be computed only *p 1 times and the 
optimization is done with respect to the following performance 
index 






(q) 2 


( 4 . 12 ) 


where 

F (q) = p(9-l) + = [f^p] 

= the resultant feedback matrix to shift eigenvalues 
of first q groups, (q=l,2* . . . ,p) 

F s the feedback matrix to shift the eigenvalues of 

Gr 

the q^* 1 group. 

The feedback matrix F & is obtained by minimizing the 
performance index (eqn. (4.12)). Thus it is to be noted 
that in this procedure the resultant feedback matrix obtained 
after shifting the eigenvalues of a group is to be optimized. 


4.9>1 Example » 

We illustrate the preceding- procedure by means of 
an example. The power system example considered in this 
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section is the same as that in Section 4.4.1. 

The performance index defined by eqn.(4.12) is to 
be minimized. The dominant eigenvalues are divided into the 
following three groups - : 

Group I : -0.014 + j0.7986 

Group II : -0.0772+ j 0.1146 

Group III : -0.0572 

By using the procedure of Section 4.9* the resultant 
optimal modal controller feedback matrix is obtained as 

”- 0.056 -0.915 1.73 1.180 -1.769 0.036 0.344 -1.725 
-0.854 -4.979 0.13 -0.058 -0.077 8.405 -1.317 -2.646_ 

The norm of the feedback matrix F as defined by 
eqn. (4.10) is 117.2. Thus this feedback matrix is better 
than the feedback matrices obtained for the suboptimal designs 
in Section 4.6.1 (Norms 233.23 and 236.07) and is comparable 
with the optimal modal controller matrix obtained in 
Section 4.4.1 (Norm,: 96.2). 

4.10 COMPARISON OF VARIOUS MODAL CONTROLLER DESIGNS 

To compare the various modal controller matrices 
obtained for the ssystems considered, a performance index 

t - I If 2 . where f . . is the element of the modal 
J ~ i=l I ij * *3 

controller feedback matrix F, is computed „ Results for the 
power system (Section 2.2.1) problem and the load frequency 
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control (Section 2.3) problem are tabulated in Table 4.1 
and 4.2 respectively. 

From the Table 4.1 it is seen that considerable 
improvement in the performance index of the feedback matrix 
is obtained by eigenvalue grouping technique of Section 4.4 
(col. (3), Table 4.1y and Col. (3), Table 4.2) as compared 
to the feedback matrix obtained without such grouping using 
method of Chapter 3 (cols.(l) and (2), Table 4.1 and 
cols, (l) and (2), Table 4.2). The subopt imal feedback 
matrix has a larger value of performance index (cols. (4) 
and (5), Table 4.1 and col. (4), Table 4.2) than the optimal 
feedback matrix (Col. (3), Table 4.1 and Col. (3), Table 4.2. 

The feedback matrix (Col. 6 , table 4.1) for the near optimal 
design (Section 4.8.1) is better than the suboptimal design 
(cols. (4) and (5), Table 4.1). The design procedure of 
Section 4.9 for large systems gives a better feedback 
matrix (col. (7), Table 4.1) then the suboptimal (Section 4.6.1) 
and near optimal (Section 4.8.1) feedback matrices (cols. (4), 
(5) and (6), Table 4.1) and the performance deterioration 
is negligible in comparison to the optimal design (col. (3), 
Table 4.1). 

4.11 OPTIMAL &ODAL CONTROL OF A LARGE ST ST EM 

The numerical examples of power systems considered 
so far, in this chapter, have systems of order less than 
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"ten and with two control inputs. The efficacy of the 
methods developed for the <f esign. of modal controllers for 
large systems can "be demonstrated effectively by consider- 
ing larger dynamic system. In this section a large system 
representing a chemical plant and described by 41 differential 
equations with 8 control inputs [73] is selected to demonstrate 
the techniques, 

4.11,1 Description of the Chemical Plant 

The following is a brief description of the chemical 
plant given in reference [73] based on the modal of Williams 
and Otto [78], The plant consists of: (l) A stirred tank 
reactor, (2) A reaction cooler heat exchanger, (3) A decanter, 
(4) A distillation column unit, which includes a reboiler 
and condenser. A brief description and data for the linearized 
model of the chemical plant considered are given in Appendix C. 
The eigenvalues of the system matrix are! 

-0,3306, -0.605, -0.815, -0.99 + 30.04, -4.905, -6.204, -6.773, 
-7.1,-9.792, -15.72, -16.2, -18.87 + jl0.48, -18.9 ± 310.78, 
-20.58, -20.6, -21.6 + 39.987, -35.19, -35.82, -41.2, 

—45.12, —45.5, —46.2, —58.59, -82.22, -82.6, —82.68, —82.71, 
-83.14, -176.9, -195.43, -350.0 (multiplicity 6), -610.25. 

It is desired to improve the system response by 
assigning the first three dominant eigenvalues [—0.3306, 

-0.605, -0.815] to the desired closed loop locations 
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C —1»5» -2*0, -2.5] by modal control feedback. 

The performance index defined in eqn.(4.12) is to 
be minimized* The modal controllers were designed to 
assign the dominant eigenvalues by applying the proaedures 
outlined earlier in Sections 3.6, 4.7 and 4.10. The three 
designs of modal controller are presented below. 

DESIGN I: Optimal Modal eohtfollor without. BigSnvr!n§ Grouping 

Applying the technique of Section 3.6, the following 
results are obtained: 

2 . = [0.30059295xl0 7 -0.31186 74 5x10 +2 -0.42868258xl0 4 

-0.42868258xl0 4 -0.96045876xl0 4 0.1l957965xl0 5 
0.14080716 0.5 2114441x1 0 4 ] ^ 

g = [-0.11025 0.10370 -0.075548 -0.07379! 0.10020 

0.11455x10“^ 0.31581xl0“ 7 0.96441xl0" 9 0.35765xl0” 9 
0.001838 0.0086633 -0.0060509 -0.0041796 

0.0051829 -0.21172xl0“ 9 0.046673 0.22110 -0.15428 

-0.10653 0.13214 0.52745xl0“ 7 1.3318 4.1223 

-1,3359 -0.44315 0.53715 0.033162. 0.18197 

-0.12381 -0.08457 0.10573 -0.017673 0.11271 
-0.079455 -0.061940 0.079465 -0.042068 0.089896 
-0.064625 -0.055312 0,072419 ] T 

The modal control feedback matrix is 
T 

F = a g 
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T,nd the value of the performance index as defined in 
oqh. (4.12) is 

J = 0.192xl0 13 


i'ESIGIT II: Optimal Modal Controller with Eigenvalue Grouping 

Using the criterion of grouping the eigenvalues 
(Section 4.2.2), the eigenvalues are divided into the 
following two groups: 

Group I : [-0,6056] 

Group II : [-0.3306, -0.8156] 


Applying algorithm given in Section 4.3 and optimizing 
':he performance index defined in eqn. (4.12) the following 


osults are obtained: 


a 



[1.0 1.0 1.0 1,01.0 1.0 -1.0 1.0] x 10 7 


g^ = [-0.64982xl0~ 3 0.66480xl0** 4 -0.21082 xlO” 4 
0.13315xl0“ 8 0.21237xl0~ 4 0.14268 xlO~ 10 

-0. 41333x10 ** 10 -0. 15117x10 "* 12 -0.32758xl0^ 11 
-0.68275xl0~ 4 -0.59891x10”*^ -0.13521xl0** 5 
0.58303xl0” 10 0.99588xl0“ 6 -0.47097xl0“ 12 
-0 • 17 41 2x10“ 2 -0.15274x10"’^ -0.34483xl0~ 4 

0 . 16911 X 10** 8 0.25398xI0~ 4 -0 . 55154xl0” 10 

- 0 . 21799 X 10* 1 -0.42013xl0“ 2 - 0 . 27945x10“ 5 
0:.67191xlO' ,8 O.lOOOlxlO” 3 -0.14058X10" 2 
- 0 . 12331 xl 0“° 3 -0.27839xl0“ 4 0 . 13482 X 10” 8 
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0.20505xl0“ 4 

-0.83610x10“^ 

-0.34459x10"' 

-0. 19144x10’ " 4 

0.10323x10“® 

0.15898x10"' 

0.63594x10“^ 

0.23987x107 5 

-0.16464x10" 

0.95468x10“^ 

0. 14816x10 “ 4 

] T 


ar d = [-2. 6024261 -6*2.4585953 -0.14609265- - 0.14609265 

-0.14606650 -0.14603529 0.14610482 1.9908551 ]xl0 7 

g^ = [0.21301xl0“ 2 -0.18491xl0" 2 0.92030xl0" 3 0.55276xl0" 3 


-0.99732x10"^ 

0.60359xl0" 12 0.10212xl0" 12 

0.36305x10" 14 

0.335lOxlO" 14 

0.345 3 2xl0" 4 

-0.13731x10"^ 

0.60353xl0" 4 

0.31309xl0" 4 

-0.48267xl0" 4 

0.73117xl0" 16 

0.88365xl0" 5 

-0.35043xl0" 2 

0.15392x10"^ 

0.79801x10"^ 

-0.12309xl0" 2 

0. 13021x1 0" 12 

0.67708x10"^ 

-0.63246X10" 1 

0.12670X10" 1 

0.33197xl0" 2 

-0.49107xl0" 2 

0 f 77986xlO" 3 

-0. 28868xl0“ 2 

0.12440xl0" 2 

0.63351x10"^ 

-0.99102xl0" 3 

0.98729x10"^ 

-0.18 5 17x10“ 2 

0.84547x10"^ 

0.46399x10“^ 

-0.75967x10"^ 

0.1l533xl0" 2 

-0.15217xl0" 2 

0.72205x10"^ 

0.41434xl0" 3 


-0.70274x10“^] 

The modal controller feedback matrix is 

]? » a (l) (g (1 V + (S (2) ) T 

and the value of the performance index as defined in 

- ? 

eqn .(4.12) is ^ 

J = 0*49xl0" 13 
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DESIGN III i Near Optimal Modal Controller With. Eigenvalue 
Grouping and Dominant Control Inputs 

The eigenvalues are divided into two groups as 
given in Section. 4.7 to identify dominant inputs, it is 
found that for Group I the dominant input is u^ while for 
Group II and Ug are the dominant inputs. In obtaining 
optimal feedback matrix the performance index defined in 
eqn.(4.12) is minimized, The following results are obtained: 

a^ T = [l.OxlO 6 0,0 0.0 0.0 0.0 0.0 0.0 0.0] 

g^ T .* [ -0 . 71913x10”^ 0. 73570x10 ” 4 -0.23331X10" 4 

0, 14735x10”® 0.23501xl0“ 4 0.15790xl0” 10 
-0.45741xl0“ 10 -0.l6730xl0“ 12 -0.36252X10” 11 
-0.75557xl0” 4 -0.66279x10""'* -0.14963xl0~ 5 
0.645 21x1 O'" 10 0.11021x10"*^ -0.52119xl0“ 12 
-0.19 26 9x10"* 2 -0.16903x10"*^ -0.38161X10" 4 
0.18714x10"® 0.28107xl0” 4 -0.6l036xl0“ 10 
-0.24124X10"* 1 -0.46493xl0“ 2 -0.30925xl0” 3 
0.74357x10“® 0.11068x10“^ -0.15557xl0“ 2 
-0.13647x10“^ -0 . 30808xl0~ 4 0.149 20xl0“ 8 
0.22692X10" 4 -0.92527x10“'* -0.38134xl0“ 4 
-0. 21186x10' “ 4 0.11423x10“® 0.17594X10” 4 
-0.70376x10“^ 0.26545x10”^ -0.18220x1c” 4 
0.1056 5x10“® 0.l6396xl0“ 4 ] 

[0.21930436x10^ 0.0 0.0 0.0 0.0 0.0 0.0 -0.86 2293x10^] 
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/ i J 

£ « [ -0.94256xl0“ 2 0.89251xl0“ 2 -0.46599xl0“ 2 

-0.30583xl0~ 2 0.5 2095x10 ” 2 -0.54128x1a” 11 
-0 . 78789X10" 11 0 . 14796xl0” 12 -0 . 81820xl0” 13 
-0.43713xl0“ 4 0.68575xl0“ 3 -0 . 31505xlQ“ 3 
-0.17322xl0“ 3 0.25500xl0“ 3 0.15664xl0” 14 
4-0.11296 x 10“ 2 0.17501x10“^ -0.80349xl0~ 2 
-0.44151xl0“ 2 0.65026xl0“ 2 0.40679xl0“ 11 
0.37896xl0“ 1 0.31894 -0.66702xl0 -1 

-0.18367x10“^ 0.26029X10” 1 -0. 12427x10 “ 2 
0.14413xl0~ 1 -0.64863xl0“ 2 -0.35050xl0“ 2 
0.52299xl0“ 2 -0.33595 x 10“ 2 0.91501xl0“ 2 
-0.43682 x10“ 2 -0. 25671x1 0“ 2 0.39953x 10“ 2 
-0.45684xl0” 2 0.74546xl0“ 2 -0.37029xl0“ 2 

-0.22924xl0“ 2 0.36864 x10“ 2 ] 

The modal controller feedback matrix is 

F = a^(g^) T + (g^) T 

and the value of the performance index as defined in 
eqn.(4.12) is 

J = 0.86xl0 13 

With either of the above modal controllers i.e. 

Design I, Design II and Design III, the locations of the first 
six dominant eigenvalues of the closed loop system are: 

-0.99 + :jl0.4, -1.5, -2.0, -2.5, -4.905 
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Davison and Chadha [73] obtained the incomplete feedback law 

u 8 = l.OxlO 7 [0.51 x 22 - 0.85 x 23 ] 
to have the locations of the first six dominant eigenvalues 
of the closed loop system as 

-0.8156, -1.0 + 0.0368, -1.6 + 37.169, -2.1826* 

The value of the performance index as defined by eqn,(4.10) 
for this feedback matrix is 
J = 0.94xl0 14 

4.11.2 Discussion of Results 

The performance indices for the Designs II and III 

are comparable. In Design III, for shifting eigenvalues 

in group I only a unique feedback matrix exists, since u^ 

is the only input. While for the second group u^ and Uq 

qualify as the dominant inputs. Therefore the optimization 

( 2 ) ( 2 ) 

is to be done only with respect to and «g elements 

of This will reduce the computational effort. 

Furthermore, in Design III, feedback is to be provided to 
inputs u-l and u Q . Therefore this design is easily implementable . 
The performance index values for the Design II and III are 
better than Design I where eigenvalues are shifted by 
applying the algorithm of Section 3*4. Considering the 
locations of the first six dominant eigenvalues of -the 
closed loop system, the performance index J for Design II 
and Design III is better than that of Davison and Chadha [73]. 
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In reference [73] > the incomplete state feedback law 
attempts to shift only one eigenvalue but there is a 

significant disturbance in the locations of several other 
eigenvalues. 

4.12 CONCLUSIONS 

In this chapter simple and computationally 
efficient algorithms have been presented to obtain optimal 
modal controllers. Grouping of dominant eigenvalues has 
been shown to yield improved results, which could be 
significant in many cases. A criterion for grouping of 
dominant eigenvalues has been given. A basis for identifying 
dominant control inputs to shift eigenvalues of each group 
is suggested. Arranging feedback to only the dominant 
inputs simplifies the control configuration and also reduces 
computational effort. Sub-optimal and near-optimal modal 
control design procedures have been suggested for design- 
ing controllers for large systems. Keeping in view the need 
for obtaining procedures to determine optimal modal controll- 
ers with reduced computational effort, a new performance index 
has been suggested. Numerical examples of modal controller 
designs for power systems have been considered in detail 
and the computational results obtained by applying the 
proced ur es developed have been discussed. To demonstrate 
the suitability and effectiveness of the procedures 
suggested in this chapter, a numerical example of the 
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design of controllers for a large system (forty first 
order system with, eight control inputs) representing a 
chemical plant has been also solved. The techniques used 
to design modal controllers for the above problem can be 
applied to power systems which are more complex than the 
one considered here. 
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CHAPTER 5 

WIDE RANGE MODAL CONTROLLERS POR POWER SYSTEMS 
5.1 INTRODUCTION 

The controller schemes discussed in earlier chapters 
were designed for the linear models obtained by linearizing 
the exact nonlinear equations describing the power system 
around an operating point. The operating conditions, 
however, change with the load demand on the system. The 
controllers designed for a particular operating condition 
may be inappropriate for other operating conditions. There 
is a need, therefore, to design controllers which give 
satisfactory system performance over a range of operating 
conditions. 

Moussa and Yu [26] considered the above problem and 
suggested incorporation of correcting feedback, in addition 
to the conventional feedback controller designed for a 
particular operating point, to compensate for small variation 
in the operating conditions^ Application of optimal linear 
regulator theory incorporating the trajectory sensitivity 
aspect has also been proposed by Elmetwally and Rao [32,33]. 
Habibullah a nd Yu [24] have suggested an iterative algorithm 
based on application of optimal linear regulator theory and 
eigenvalue search technique* 

In this chapter a scheme to determine feedback, in 
addition to the modal controller designed for a particular 
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operating point has heen developed. This additional 
feedback; obtained from eigenvalue sensitivity anal ysis 
reduces the variations in the closed loop system poles with 
change in operating conditions. This feedback is a function 
of the system state variables and those parameter of the 
system that are subjected to change. This technique has the 
advantage and simplicity of the modal control approach; 
furthermore it enables realization of the desired system 
performance over a range of operating conditions. It has 
been applied to a power system problem and the computational 
results have been given. The problem of implementation of 
this scheme has also been discussed. 

The proposed controller may be called an Adaptive 
Modal Controller or a Wide Range Modal Controller. 

5.2 DESIGN OF WIDE RANGE MODAL CONTROLLER 

5.2,1 Conventional Modal Controller j 

th * 

Consider the linear, time-invariant, n order system 

with m inputs j 

x = Aq, x + B 0 u (5.1) 

where the suffix 'o' denotes the value at the nominal [ 

i 

operating point 

Jt has been shown in Section 4 ♦ 2 , that the dynamic 

response of the nominal system can be improved by the modal | 

I 

(5.2) | 


control feedback law, 
u = F 0 x, 



33 

(5.3) 


F o * E d^(g 0 ^) T 

< 1=1 

which shifts the 'V dominent eignevalues of A 0 , by state 
feedback to the desired closed loop locations without 
disturbing the locations of the remaining eigenvalues. It 
is assumed that the dominant eigenvalues are divided into 
*p* groups. 

5.2.2 Wide Range Modal Controller 

Consider now the perturbed system 

x as (A 0 + AA) x + (B q + A B) u (5.4) 

where A and B are variations in the system and control 
matrices due to the perturbations in system operating conditions. 

The feedback law 

u = Up = P 0 2 (5.5) 

is such that for the unperturbed system, the *A* dominant 
eigenvalues of the system are assigned to the closed loop 

locations, P 1 ,P 2 > » P fc » without disturbing the locations 

of the non-dominant eigenvalues. However, perturbations i# 
system parameters will, in general, result in perturbations 
of the above closed loop eigenvalues. The objective in the 
design of a wide range modal controller is to provide 
additional feedback such that the first order variations of 
the above eigenvalues with respect to system parameter 
variations is rendered zero. For this purpose the eigenvalue 
sensitivity reduction technique given by Tzafestas and 
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Paraskevopoulas [79] for single input systems is modified 
to make it applicable for multi— input systems. 


The dominant eigenvalues are divided into p groups 
for the eigenvalue assignment problem. How, consider a 
feedback, in addition to the given in eqn. (5.5), of the form 

A u = ( E ) A6 T H x (5.6) 

1=1 ~ 

where 0 is the vector of those elements of system matrix A 0 

and the matrices B Q cS^t B 0 a^, , B Q which undergo 

perturbation, and A© is the resulting perturbation of the © 
vector. Let © be of dimension r. Then H is an r x n matrix 
and the i^* 1 row of H matrix gives the correction required to 
reduce the variation in closed loop pole locations due to 
perturbation in the i th element of ©. The elements of H 
matrix are derived as follows . 


Let the eigenvalues of the closed loop disturbed 
system be denoted as X ^ (^©) , (i=l,2, . . . ,n) . Expanding 
x 9 (©) in a Taylor series expansion around nomi n al 
value (X J (®)] 0 

3 A ?(©) 3A?(©) 

li <®>3o + [ ’ 


*i(S) 


Tsr. 


i ♦ * * 9 


« |{e) 

^6 ]*S 


+ high order terms (5.7) 

when the correction given by eqn. (5.6) is added to the 
control law obtained from eqn. (5 &S , the resultant closed 

„ Z 

loop system matrix becomes 



Then 


A ° L " A + q=l ~ (<l) + a - T ••• 


(5.8) 


It has been shown in Appendix D, that the eigenvalue 
sensitivity, [80], for the above system is given by 


3X C 
i 


-1 


<atfr->o = C T r R °<»i>]o~ [ R °(>i)] 0 * [(Ae .) 0 + \ (4 q) ) 0 

3 Q. 3- 3 

(g 0 (q) ) T + b 0 h* + H.O.T.] , (i=l,2,...,n), 


( 3 — 1> 2, . . . ,r) 


(5.9) 


where 


b a | * u) 

q=l “° 


^o 


(q) ^ B rt a (q) ; 


Aq and bQ 
D 3 


(q) 


,(q) 


with 


are the partial derivatives of A and b 

d • - d 

respect to element ©^ of © $ R c (X^) = Adjoint [IX^ - Aq - 

7; b (g_^^)^] ; h* is the j^* 1 row of H and the suffix 
q=l ~° - ” J 

’o’ denotes that the expressions are evaluated at the nominal 

J3 

operating conditions. X? , (i = l,2,...n) denotes the desired 
closed loop pole locations corresponding to ^ g, * 

X S,+ !>••"•» X n respectively and the operator & is as defined 
in Appendix D. 
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riiL 

r\ 


The second term of eqn* (5.7) is zero if 


■ ) 0 ss 0, (j = l,2,...,r) and (i=l,2, . . . ,n) . Then, from 


eqn. (5.9) 


[R°(lf)] 0 9 [b 0 hj 1 ] = - tR C (J ± )] 0 * CU 6 .) 0 + 


? (£o. (9) ) 0 (& (4) >^ 

q=l 3 


rCR c Cxf)3 0 b 0 ] T h 3 = - [R°(x ± )] 0 9 [(Ag ) 0 + 


? (b 9 (q) ) 0 (g 0 ( ' l) )' r ] (5.10) 

q=l 3 


Q h^ = — e j t ( j=l» 2> t • • »n) 


(5.11) 



n-vector 
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The solution for hu , (j=l,2, .. .,r) is obtained 
by solving the equations (5.11)* If Q is singular, the 
equations in the set (5.11) corresponding to the rows of 
Q which are linearly independent are retained and the 
remaining equations are discarded. A solution for hj can 
then be obtained. H 0 wever the solution obtained is not 
unique . 


5.3 WIDE RANGE MODAL CONTROLLER POR A POWER SYSTEM 
5.3.1 System Model 

The power system considered is the same as that 
in Chapter 3 Section 3.6.1. The numerical data for this 
system is given in Section 2.2.1. The desired closed loop 
pole locations for the nominal system are the same as those 
given in Sections 3.6.1 and 4.4.1. 

It is desired to design a wide range modal controller 
so that the deviation in the closed loop eigenvalue locations 
for a wide range of operating conditions is small. The 
operating conditions are given by the real power P and the 
reactive power Cl delivered by the generator. The elements 
of the system matrix A, which depend on P and Q are 

a 2l’ a 23* a 31 8,41 (See Secti ° n 2 ' 2 * 1) * The elements 


CJL ty y — . 

of the control matrix B are independent of the operating 


conditions. 
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5.3.2 Results 

Two modal controller designs are considered for 
shifting the five dominant eigenvalues of the uncontrolled 
system to the required closed loop locations at the nominal 
operating conditions. In design I all the five dom inant 
eigenvalues are considered to be in one group and the pole 
shifting procedure given in Section 5.2.1 is applied. In 
design II, domin^it eigenvalues are divided into two groups. 
The first three domii^nt eigenvalues are in group 1 and the 
next two are in group 2. The pole assignment procedure is 
then applied to the two groups ag indicated in Section 
5.2.1. Using sensitivity considerations and the procedure 
of Section 5.2.2. additional feedback controller is designed 
for each of the above cases. Closed loop system pole 
locations for the above two designs are computed for various 
operating conditions. The results are tabulated in Tables 
5.1 to 5.5. The operating conditions (real and reactive 
power delivered by the machine) considered also include the 
cases considered by Habibullah and Yu, [24]. In table 
5.5 closed loop system eigenvalues for the controllers 
with and without additional feedback for certain operating 
conditions are compared. 

DESIGN I 

Modal Controller Design [see Section 3.6.1] 
a = [-25.1027 , 94.6869] 1 
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£o = [°.055» 0.104, -0.165, 0.267, -2.528, 8.3466, -0.65, 

0.935] 1 

b Q = [0.0, 0.0, 0.0, 0.0, -25.1027, 0.0, 94.6869, 0.0] T 


Design of Additional Feedback 

A 0 = [A a 2]_ » A a 23» A a ^i » 4a^-^] 

0.104 0.640 -0.314 0.201 -2.170 10.130 -0.575 2.216 

0.164 0.026 0.677 0.593 -3.438 10.661 -0.917 0.760 

-0.097 -0.315 0.460 0.201 0.328 1.291 0.087 1.045 

-0.412 0.201 0.201 -0.532 4.992 -5.888 1.323 3.611 

The wide range modal controller feedback matrix is 
P wr = So §o T + \ ^ H (5 ’ 12) 


DESIGN II 

Modal Controller Design [see Section 4.8.1] 


cf 

H 

II 

[ 0.0 i.of 

0.183, 8.38, -1.32, -2.678] 1 

ii 

^ — y 

H 

o 

t*0l 

[0.881, -5.028, 0.086, -0.1, 

£o (1) = 

[0.0, 0.0, 0.0, 0.0, 3.122, 

0.0, 3.296, O.O] 1 

= 

[1.0 0.0] T 


go (2) - 

[0.481, 0.184, 1.818, 1.204, 

-1.766, -6.566, 0.093, 



-2.774] 1 

II 

C\J 

K— * 

c 

[0.0, 0.0, 0.0, 0.0, 0.3749, 

0.0, 3.22* O.O] 1 

&© . •*. ~ 

o.o.,. o.o, o.o, 1.0, o 

.0, 1.0, 0.0 ] 

DesJ.gn 

o-f Additional Feedback' 

; * v . 


A© *= [A 

a 21» Aa 23 f Aa 31* A 41 J 
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H 


-9.039 -5.137 -1.474 -0.339 0.296 29.194 -0.296 12.670 
0.135 -1.319 -0.546 -0.138 0.120 -6.650 -0.120 -0.39 6 
-1.019 -1.474 -2v'69l -0.526 0.400 14.257 -0.400 4.981 
-0.269 -0.339 -0.526 -0.122 -0.099 3.174 -0.099 1.060 

Tlie wide range modal controller feedback matrix is 





+ b A9 T H 

—A — 


( 5 . 13 ) 


5.3.3 Discussion of Results 

The results given in Tables 5.1 to 5.5 show that the 
variations in the closed loop pole locations have been 
rendered small over a fairly wide range of system operating 
condition by the Wide Range Modal Controller scheme. 


Examination of the two designs indicates an interesting 
aspect of the additional feedback in its relation to the 
nominal modal control feedback. The correction required to 
the nominal control could be viewed as being indicated by the 
magnitudes of the elements of the matrix H. larger the 
magnitudes of the elements of the matrix H, larger the 
correction provided by the additional feedback. It is 
observed that the correction provided by additional feedback 
is more in Design II than that of Design I. This can be 
explained as follows: The nominal modal control feedback 
gains in Design I are larger than those In Design II. 

Therefore the relative variations in the elements of the 



101 


TABLE 5.1 : CLOSED LOOP SYSTEM EIGENVALUES FOR 
DIFFERENT REACTIVE LOADINGS , 


S.No . 

Q (p.u) 

Closed loop system eigenvalues with 
Via© ..Range Modal '.Controller... _ 

1 

0.7 

-13.71, -1.553, -0.7122 + jl.272, 



-0.1814 + j0.4802, -0.2858, -0.1559 

2 

0.6 

-13.71, -1.4, -0.7028 + jl.176. 


-0.2637 + jO.5145, -0.2834, -0.1655 

3 

0.5 

-13.71, -0.6875 + jl.072. -1.222, 
-0.3646 + 00.5589, -0.2809., -0.1743 

4 

0.4 

-13.71, -0.6554 + jO.9384, -1.001, 
-0.5046 + 0*0.6443, -0.2783, -0.1820 

5 

0.3 

-13.71, -0.7979 ' + 0*0.7151, -0.7713, 
-0.4751 ± 0*0.871 5, -0.2761, -0.1883 

6 

0.2 

-13.71, -0.9092 + 00.7171, -0.6508, 
-0.4226 + jO.9025, -0.2745, -0.1927 

7 

0.1 

-13,71, -0.9481 + jO.7245, -0,6072, 
-0.4048 + j0.9137, -0.2736, -0.1950, 

8* 

0.034 

-13.71, -0.9539 ± 0*0.7255, -0.6006, 
-0.4021+J0.9155, -0.2736, -0.1954 

9 

0.0 

-13.71, -0.9523 + o0.7251 r -0.6027, 
-0.4028+ 0*0.9151, -0.2735, -0.1953 

10 

-0.1 

-13.71, -0.9332 + j0.7212, -0.6244, 
-0.4114 ± 00.9094, -0.2739, -0.1940 

11 

-0.2 

-13.71, -0.8924 + jO.7167, -0.6733, 
-0.4290 + d 0.8971 +0.2739, -0.1915, 

12 

-0.3 

-13.71, -0.8263 + 0.7184, -0.6944, 
-0.4575 ± 00.8768, -0.2737, -0.1884 

13 

-0.4 

-13,71, -0.726 + o 0.7478, -0.8584, 
-0.5068 + 0*0.8379, -0.2731, -0.1848 
-13.71, -0.6622 + o0.8922, -0.9746, 

14 

-0.5 

-0.6622 + 00.8922, -0.2721, -0.1808 


Nominal operating condition P Q =0.735, Qq= 0.034, v t -1.05 
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TABLE 5.2 : CLOSED LOOP SYSTEM EIGENVALUES POE 

DIFFERENT REACTIVE LOADINGS DESIGN II 


S.No 


1 

2 

3 

4 

5 

6 
T 
8 
9 

10 

12 

*3 

14 




• Q ( P • ) 


Closed loop system eigenvalues with 
Wide Range Modal Controller 


0.7 


0.5 


0.4 


0.3 


0.2 


0.1 


0.034 


0.0 


- 0.1 


- 0.2 


-0.3 


—0.4 


-0.5 

-0.7 


-13.72, -0.9737 + d0.7154, -0.6231, 
-0.3966 ± jO.9389, -0.3013, -0.1175 . 
-13.71, -0.9672 + jO.7154, -0.6033, 
-0,5988 + 00.9275, -0.2902.,. -0.1587 . 
-15.71, -0.9649 + 00.7159, -0.5979, 
-0.3995 ± 00.9235, -0.2846, -0.1739 . 
-13.71, -0.9630. ± o'0.7167, -0.5953, 
-0.3999 + 00.9204, -0.2797, -0.1851 . 
-13.70, -0.9616 + JO. 7177, -0.5954, 
-0.4001 ± 00.9181, -0.2761, -0.1922 . 
-13.70, -0.9605 + dO.7191, -0.5978, 
-0.4000 ± JO. 9163, -0.2743, -0.1951 . 
-13.70, -0.9600 + JO, 7201, -0.6004, 
-0.3999 ± J0.9154, -0.2740, -0.1951 . 
-13.70, -0.9597 ± o'0.7207, -0.6021, 
...-0,3998 + jO.9150, -0.2741, -0.1945 . 
-13.70, -0.9592 + jO.7225, -0.6080, - 

-0.3993 ± d 0.8140, -0.2748, -0.1714 . 
-13.69, -0.9588 + d0.7226, -0.6154, 
-0.3982 + d 0.9132, -0.276, -0.1866 . 

-13.69, -0.9581 + d 0,7 209, -0.6239, 
-0.3979 ± d0.9125, -0.2771, -0.1810 . 
-13.68, -0.9586 + d°*7239, -0.6335, 
-0.3970 +jO,9121, -0.2777, -0.1748 . 

-13.68, -0.9589 + d0.7320, -0,6442, 
-0.3960 + d 0.9116, -0. 2778, _ -0.1683 
-13.67, -0.9604 + d0.7376, -0.6691, 

-0.3937+ 30*9114, -0.2758, -0.1541 


Normal operating condition Pq-0,735, Qq-0.034, 


.05. 
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S.No, 


TABLE 5.3 '• CLOSED LOOP SYSTEM EIGENVALUES FOR 

DIFFERENT OPERATING CONDITIONS .DESIGN I 

• Q Closed loop system eigenvalues with, 

p.u, p.u. Wide Bange Modal controller 


-0.3168* -0,1811 + 3 0.4649, -0.1746 . 

• 1 ^ 20 0 ^ 4 -13.71,-1.2141,-0.6787 + 31.103, 

-0.3405 + 30.5036, -0.3150, -0.1876 . 

!• , ,, n oc *13.71, -0.8307, -0.6283 ± 30.5992, 

X * It) U • <£j 

-0.5945 ± 30.9259, -0.3078, -0.1975 . 

r o 952 o 015 ~ 13 * 71 ' “°-9974 + 30,7323, -0^5357, 

-0.3875. ± 30.9312, -0.2731, -0.2028 . 

► 0 7 -0.15 “ a3 - 71 » “0*9514 ± 30.7245, -0.6026, 

“ S -0.4035 + 30.9149, -0.2739, -0.1953 . 

' 0 5 -0.225 - 13 * 71 ' - 0 *9658 + 30.7192, -0.5631, 

* -0.4037. + 30.9226, -0.2804, -0.1991 . 

► 03 -0 256 - 13 * 71 ' - 1 - 055 ± 30.7414, -0.4153, 

-0.3696 + 30.9539, -0.3081, -0.2097 . 

n . 0 02‘56- 13 - 71 ’ _:l - 050 ± 3°- 7436 ’ -°' 4618 * 

° ' -0.3691 + 30. 9518, -0.2710, -0.2115 . 

. q „ n lt - -13.71, -1,001 ± 30.7336, -0.5332, 

U * y7 -0.3855 ± 30.9324, -0.2717, -0.2032 . 

-13.71, -1.112 + 30.7712, -0.4064, 

1.15 *0.25 _o.3404._+ 30.9755, -0.2355 + 30.0419. 

-13.71, -1.143 ± 30.7852, -0.3845, , 

1.20 -0.34 -0.3262+ 30.9891, -0.2290 + 30.06 . 

-13.71, -1.176 + 3 0.8014, -0.3710, 

1.25 -0.45 _o.3117 + 31.004, -0,2174..+ 30.0754 . 

-13.71, -1.307 + 30.8786, -0.3339, 

1.4 -0.6 -0.24 + 31*063, — 0.1771 + 30.1121 • 

-13.71, -1.382 + 30.9298, -0.3237, 

1.5 -0.65 -0.1935 + 31.H, -0.1528 + 30.1252 

operating conditions considered by Habibullah 

and Yu [24] 
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TABLE 5.4 : CLOSED LOOP SYSTEM EIGENVALUES EOR 


DIFFERENT OPERATING CONDITIONS DESIGN II 

I 


S.No. 


1 

p.u. p.u. 

1.5 0.65 

Wide Range Modal Controller 

-13.67, -0.9782 + 30.7073, -0.7053, 

2* 

1.25 0.45 

-0.4069, + 3 0.9078, -0.3439, -0.0526 . 
-13.68, -0.9720 + 3*0.7049, -0.6474, 

3* 

1.20 0.34 

-0.4025 + o' 0.9046, -0.3165,. -0.1332 . 
-13.68, -0.9702 ± jO.7051, -0.637, 

4* 

1.15 0.25 

-0.4017 ± 3 0.9031, -0.3048, -0.1562, 
-13.68, -0.9692 + 30.7062, -0.6318, - 

5 * 

0.95 0.015 

-0.4009 ± j0.9033, -0.2942, -0.1712,. 
-13.68, -0.9653 + j0.7144, -0.6246, 

6* 

0.7 -0.15 

-0.3989 ± 3*0.9086, -0.2739, -0.1900 . 
-13.7, -0.9588 + j0.7217, -0.5986,.. 

7 * 

0.5 -0.225 

-0.4 + .3*0.9159, -0.2742, -0.1952 . 
-13.71, -0.9522 + 3*0.7279, -0.5804, 
-0.4018 + 3*0.9198, -0.2769, -0.1975 . 

8 * 

0.3 -0.256 

-13.73, -0.9491 ± 3 0.7269, -0.5479, 
-0.4053 ± 30.9289, -0.2675, -0.2131 . 

9 

0.7 0.15 

-13.70, -0.959 + 3*0.7211, -0*5971 , 
-0.4001, + 30.9163, -0.2739, -0.1958. 

10 

0.95 _0.15 

-13.68, -0.965 ± 30.7148, -0.6266, 
-0.3987 + j 0.9085, -0.2732, -0.1896 . 

11 

1.25 -0.45 

-13.65, -0.9806 + j0.7031, -0.6934, 
-0.3961 . + JO. 09069, -0.2568, -0.1671 . 

12 

1.4 -0.6 

-13.63, -0.988+ 3 0.6907, -0.7141, 

-0.3974 + J09055, -0.2488 x -0.1648 . 

13 

1.5 -0.65 

-13.63, -0.9895 ± jp.6818 -0.718,: 
-0.3997 + 30*9005, -0.2476, -0.1707. 


operating uv 

and Yu [24] 
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closed loop system matrix due to the change in operating 
conditions are smaller in Eesign I. Hence the amount of 
correction required in Design I is smaller. 

5.4 COMMENT ON PHYSICAL REALIZATION 

In the implementation of the wide range modal 
controller, state feedback proportional to the c han ge in 
system parameters is to be..- provided. This can be achieved 
in the case of the power system considered in Section 5.3 
using advanced techniques in adaptive control such as model 
reference adaptive control schemes which can be implemented 
on line. Por the problem of power system considered, the 
real and reactive power delivered by the machine are easily 
measurable quantities. 

5.5 CONCLUSIONS 

A simple and effective method of designing wide range 
modal controllers for the power system is presented in this 
chapter. The controller so designed reduces the variations 
in the closed loop system eigenvalues to variations in 
operating conditions over a fairly wide range. This assures 
satisfactory system performance over a wide range of system 
operating conditions. It is observed that if the feedback 
gains in the modal controller design for nominal operating 
condition are reduced significantly by optimization, then 
large correction to the nominal control is required when the 
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operating condition changes. These observations are in- 
agreement with the result in classical control theory, 
which, states that larger the negative feedback, tha more 
insensitive is the system performance to variations in the 
system parameters [ 81 ], 



CHAPTER 6i 

ACCESSIBLE STATE FEEDBACK MODAL CONTROLLERS 


6.1 INTRODUCTION 

The modal controllers designed in earlier chapters 
require feedback of all the state variables of the system. 
This may be difficult to achieve in practice. If the system 


is observable, it is possible to reconstruct the state 
variables to the required degree of accuracy from the 
system outputs with the help of 'observers’ [82]. The 
reconstructed state variables can then be fed back to the 
conventional feedback controllers. However, the disadvantage 
of this technique, especially in high order systems, is the 
additional dynamics due to the observers. Davison [83] 
obtained results regarding assignment of *( <n) open loop 
poles to desired closed loop locations by linear output 
feedback. In this procedure, the unassigned poles may get 
disturbed. This may cause the closed loop response to 
unsatisfactory. Fallside and Seraji [84, 85] proposed a 
method for the output feedback pole assignment problem with 
unity rank matrices. In this case also the locatio 

unassigned poles may be disturbed. 

,, nf modal controllers 

In this chapter the desig 

. . nut rut feedback to assign the 

utilizing accessible state or P 

, +n desired locations, without disturbing the 
dominant poles to desirea 
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locations of the remaining poles of the open loop system, 
has been considered. The possibility of optimizing feedback 
gains has been investigated, numerical examples have been 
given to illustrate the procedures developed. The cha pter 
also includes design of a dynamical modal controller scheme. 
The schemes developed have been applied to a power system 
problem. 


6.2 MODAL CONTROLLERS USING ACCESSIBLE STATE FEEDBACK 

Consider the linear, n order system with m control 

inputs described by 

x = A x + B u (6.1) 

Let 

where 

x-^ = the r-vector of inaccessible states. 

It is required to find a feedback law u = F 2 x 2 to 
shift the eigenvalues ^2 , *** , \ P reassi § ne<a - locations 
p l» p 2 , . . • » P respectively without disturbing the locations 
of the remaining eigenvalues of the uncontrolled system. It 
is assumed that the modes ^ to X^are controllable. 


Let F 2 
F 2 = 2, 


be constrained to be of. the dyadic form: 

T (6.2) 

2.2 


where a is an m-vector and the feedback gain vector ^ is 


of dimension (n-r). 



Ill 


The gain vector g of dimension n for the complete 
state feedback pole assignment problem is given by eqn.(3.6). 

£ = V l2. (6.3) 

where 

- th 

I component of the ^“Vector c 

= [V<3Ci T B a)] (6.4) 

a i 

end kj » [ ^ (p ± - Xj)]/[ n (X^X^)], ( j=l,2,. .. ,5. ) 

(6.5) 


It is to be noted that for the pole assignment problem, 
CjL 5^ (i=l,2 f .*♦, Jl ) 

let 

v L . [ i* , v 2 t ]* 

§ T ” [ Si ' Kg 1 ' 3 


and 




where "V-^, g_^ and c-^ are of dimension rx Z, r and s respectively, 
s being the rank of V 1 , If the state vector x ± is not to be 
fedback, then g ± = 0. Therefore, from eqn. (6.3) 

V 1 c = 0 (6.6) 

A necessary and sufficient condition for non-trivial solution 
c ifceqn. (6.6) is s <fl, -1. This implies that r<A -1, that 
is, the number r of inaccessible states should be less than 
the number of eigenvalues to be shifted. Since 
rah% [V 1 ] » s, ( A-s) elements of c can be arbitrarily selected 
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and the rest of the elements determined from eqns. (6.6). 

It is required that none of c.'s determined should be zero. 
In many problems, solution for c satisfying the above 
conditions can be obtained. 

Equations (6.3) are rewritten below for obtaining a: 

V L B St = £ (6.7) 

th 

whore the i element of the & -dimensional z vector is 

z i ~ ( k i/°i) # Let rank [V = rank [Y^ B ! z], Then 
a solution for a can be obtained. Then F 2 is obtained after 
determining g from eqns. (6.2). 

It is to be noted that if s < r and/or &< m, the 
freedom available in tthe choice of c and/or a can be 
utilized to obtain optimal feedback gains. Furthermore, if 
it is found that the solution of eqns. (6.7) does not exist 
(this is more likely to happen if £ is large), then the 
eigenvalues to be moved may be divided into groups and the 
procedure may be applied repetitively to the groups. 

6.3 EXAMPLES 

Three numerical examples are considered in this 
section. The first two are simple problems to demonstrate 
the suitability of the procedure developed in Section 6.2, 
for designing accessible state feedback modal controllers. 

In example 3» accessible feedback modal controller design 
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for a power system is given* 

6,3.1 Example 1 : Accessible Feedback Modal Controller 

for a Simple System 

The system considered has the following system and 
control matrices: 



0 10 


0 0 

A = 

-6 -5 1 

i BZ== 

■ 0 1 


0 0-1 


1 0 




r— 


The measurable states are 


x 2 = [ x x , x 2 ] T 

The eigenvalues of uncontrolled system are -2,-3 and -1. 

A modal controller which utilizes feedback from the measurable 
states x 2 is to be designed to shift the open loop eigenvalues 
-2 and -3 to locations -4 and -5. 

Solution: 

r*=i; i= 2 ; i = 2 ; n = 3 ; 
k 1= -6 J k 2 = 2 i c i = for c 2 = 1.0 
a * [2.5, 1.0] T i g = C'r4. -4 » 0] T 

and the modal controller feedback matrix is 


-10 -10 
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6.3.2 Example 2 • Optimal Modal Controller with Accessible 

State Feedback 

Consider the example 1, in which the B matrix is 


changed to 






0 

0 

1 

B = 

0 

1 

1 


1 

0 

0 

Solution: 





Since &< m, the elements of a vector are selected to 
minimize the performance index 

3 2 ^ 2 

J “ ill a i + 3-1 


Then, 

g. « [-1.915, -3.566, 1.0 f 

g * [2.56, 2.56, 0.0] T 

and the modal controller matrix with measurable state 
feedback is 


-4,9 

-4.9 

-9,13 

-9.13 

2.56 

2.56 


6.3.3 Example 3s Accessible State Feedback Modal Controller 

for a Power System 

Ihe power system considered is the same as that given 
in Sections . 232.1 and 3.6.1. 


The complete state vector is ■ 

X- lit, 4s ' iS f ,ih] 
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The state variable represents the flux linkages of the 
field winding. This quantity cannot be directly measured. 

A modal controller which does not require feedback of A 1 ? ^ 
is to be designed to shift the dominant eigenvalues 
[-0.0114 ± 30.7986, -0.0572, -0.0772 + $0.1146] to the 
closed loop locations [-0.4 + 30.915, -0.6, -0.96 + 30.72]. 
Since r — 1, we take £ =2. The number of dominant modes is 
five. In order to have 3 groups with two eigenvalues in 
each group, one of the non-dominant eigenvalues is also 
shifted, but only very slightly. Thus the eigenvalues to 
be shifted are divided into following three groups: 

Group I * [-0.0114 ± j 0.7986] 

Group II 4 [-0.0772 ± 30.1146] and 
Group III : [-0.0572, -0.1953] 

and the respective desired closed loop locations are, 

[-0.4 ± 30.915], [-0.96 + 0.72] and [-0.6, -0.196]. 

The procedure of Section 6.2 is applied in a repetitive 
manner* The resultant modal control law is found to be 

. • ~ 12.3452 -461.1443 0.0 -0.1196 -10.6961 -3801.6890 

F is® ' 

-40,6185 -20.9747 0.0 -0.0960 -0.6996 -222.0193 

-124,4492 -1295.2454 

7.6098 -73.5784 

Note that the column corresponding to A* f in the above 
modal controller feedback matrix has zero entries. 


(6.9) 



Torque angle in radians AS 



FlGtTHE 6.1 SYSTEM RESPONSE FOR AS 


Torque angle in radiance A6 
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It may be possible to have a better feedback matrix 
in term of the performance index given by eqn. (3«8), if 
instead of the non-dominant eigenvalue -0.1953, some other 
eigenvalue from the non-dominant set is shifted. For 
example, when the eigenvalue -13.4 is shifted to -12*33 
alongwith the dominant eigenvalues, the feedback matrix is 


110.5713 1282.31 0*0 16.7521 -9.7971 -312.3934 
0.9469 13.99 0.0 0.2001 -0.1447 -4*6833 

-131.6718 218.8719~ 
8.0858 2.7420 


( 6 . 10 ) 


The feedback gains can be further reduced if eigen values at 
location -13*4 is shifted towards the right half plane without 
appreciably affecting the closed loop performance. When it 
is assigned -5*28, the feedback matrix is 

~46.9245 528.2866 0.0 6.7510 -4.060 -145.8445 

F lies 

0.4722 8. 3518 0.0 0.1255 -0.102 -3*4411 

43*6062 85.8360] 

( 6 * 11 ) 

9.3930 1.7499_ 

The closed loop system responses with the above 
controllers and the initial conditions considered in 
Section 3-6.1 are given in Figure 6.3. 


6.4 OUTPUT FEEDBACK MODAL CONTROLLER 

Consider the system eqn. (6.1). The p outputs of 

the system are given by the output equation 
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Lot rank [D] = p. It is required to find a feedback law 

a = p 0 y (6.13) 

to shift X 1 , X 2 ,..., X^ to preassigned locations P- L ,P 2 , .. . ,p £ 
without disturbing the locations of the remaining eigenvalues. 
It is assumed that the modes X. toX are output controllable. 

Xj 

Let F Q be constrained to be of dyadic form 

E o “ (6.14) 

whore the gain vector p is of dimension p. 

Now, consider the case when it is possible to generate 
g of eqn. (6.3) » in the row space of D, that is, 

g s= D® £ (6.15) 

Then 

a » a (D T p) T x (6.16) 

The control law (6.16) utilizes only output feedback. From 
(6.3) and (6,15) we have 

jp p as c (6,17) 


or 



css 0 


(6.18) 


where D » [D T : - V L J is an nx(p+A) matrix. Let rank [5] * r. 
Then a necessary and sufficient condition for non-trivial 
solution for [P T s c T ] in eqn. (6.18) is 


(6.19) 


p + % > r 
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It is required that none of the o 1 's determined 
should be zero. In many problems, solution for o satisfying 
the above criterion X'found. Equations (6.7) can now be 
solved to obtain a . The output feedback nodal controller 
matrix F Q can then be obtained from eqn. (6. 13). 

T 

If rank B] < m, then the freedom available in 
the choice of gg may be utilized to set some of the elements 
of gc_ to bo zero. This implies that feedback is arranged 
to only some of the available inputs, resulting in a 
simpler control configuration. 

6.4.1 Optimal Modal Controller 

If p+ & - r > 1 in eqn. (6.19) and/or rank [Vjj: B] < m, 
then freedom available in the choice of c and/or a_ in 
eqn. ( 6,18) and eqn.( 6.7) can be utilized to obtain optimal 
feedback gains 1 The following performance index may be 
selected for feedback gain minimization.* 

J = E E (6.20) 

i=l j=l J 

6.5-. ALGORITHMS 

6.5.1 Algorithm To Design Output Feedback Modal Controller 
£ ) Determine the reciprocal eigenvectors of the system 

matrix A, corresponding to 1 & ’ eigenvalues 
( ) to be assigned. Also check that 

modes X to X are output controllable, 

x j, 
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Jli) Calculate k^, (j=l,2,,.., SL ) from eqn. (6.7). 

lU) Compute the rank of the matrix D (eqn. (6.18) ). If 
condition eqn. (6. 19) is satisfied, then go to step 
(iv), otherwise stop. 

iv) Obtain the solutions for the vector (3 and c from 
eqns. (6.18) such that c i ^ 0, (i=l,2,..., $,). 
r) Compute [k.j/cJ, (i=l,2,..., z ). 

vl) Jj ind the values for (i=l,2,...,m). by solving 
eqns. (6,7). 

vii) Obtain the modal controller matrix F from eqn.(6.14). 

viil) Stop. 

J.5.2 Algorithm to Design Output Feedback Optimal 
Modal Controller 

i) Determine the reciprocal eigenvectors of the system 
matrix A, corresponding to 1 i 1 eigenvalues 
( \ lt X , X > to be assigned. Also check that 

X, to X are output controllable, 
j. ^ __ 

ii) Compute r the rank of matrix [D] eqn. (6.18). If 
p+ A -r eqn. (6.19) is not greater zero, then stop. 

Let p+ a -r = r^. 

iii) Calculate k,., (3=1,2,,. t , i ) from eqn. (6.15). 
iv) Set iteration count n=l. 
v) Set objective function value J - 0, where the 
objective function J is of the form (6.20) and 
superscript n denotes the iteration count. 
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vi) Choose suitable values for r ± elements of the c 
and let these elements constitute the vector Cg. 

Obtain the solution for £ and the remaining element 
of c from eqns. (6.18). The vector c. 2 should be such 

that £ 0, (i=l,2, . . . , i ). If the above requirement 
cannot be met then stop. 

vii) Compute z ± = [k^], (i=l,2,..., i ). 
viii) Find r 2 , the rank of matrix [v£ B] in eqns. (6.7). 

If rank [vijj? B] ^ rank [V^ B : z], modify c 2 in step 
^i) *to make it equal. If this is not possible then 

stop. 

ix) Choose suitable values for (m-r 2 ) elements of a. 

Lot these elements form a vector a 2 . Solve eqns. (6.17) 
for the remaining elements of a. 

x) Calculate the new objective function value j( n+1 ) 
from eqn. (6.20). If I J^ n+1 ^ - I is less than 
a specified small positive number, go to step fevi), . 

xi) Change £ 2 and a 2 vectors as c 2 new - 9. 2 + s ( A 2^ 
and &2 now = <% 2 + s (A a 2 ) respectively, where s is 
the step length and A c, 2 and Aa 2 are the negative 
gradients of the Lagrangian function with respect, 
to c 2 and a g respectively. 

xii) Solve eqns. (6.18) for £ and the remaining elements of 
£ . Check, whether c. / 0, (i=l,2,..., * ). Otherwise 
modify appropriately the step length s in stepCxi) . 
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xiii) Repeat step vii. 

xiv) Solve oqns. (6.17) lor £ . if solutions for £ 4o 
not exist then change suitably the step length 
s in step(xi). 

xv ) Increment the iteration count n by one and go to 
st ep(x). 

xvi) Obtain the modal controller feedback matrix F 

o 

from oqn. (6.14). 
xvii) Stop, 


6.6 NUMERICAL EXAMPLES OP THE DESIGN OP OUTPUT 
FEEDBACK MODAL CONTROLLERS 

'Throe numerical examples are considered in this 
section to demonstrate the suitability of the algorithms 
developed in the previous section for designing output 
feedback modal controllers. In example^ 2, the feedback 
gains are optimized. 


6.6,1 Example 1 : Output Feedback Modal Controller 

for a Simple System 

Consider the system (6.8) of example 6.3.1. Let 

the output matrix be 


D 


1 1 0 

0 1 0.5 


The output feedback modal controller is to be 
designed for the desired closed loop, pole locations 



mentioned in example 6 . 3 . 1 . 
Solution; 

In this problem; 
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P = 2 > &= 2, m=2, r=2 and p+ £ -r = 2 > 0. Taking 

c 2 = 1 * 0 , we obtain from eqns. (6.18) 

C-J^ sa 1.0 

« [ 1 ,- 1 ] 

Then 

a T =[82] 

and the feedback matrix 


P jgg 
" 0 


8 

2 


-8 

-2 


6.6.2 Example 2 : Optimization Problem 


Consider the example 6.6.1, except that B is 


changed to 


B = 


0 0 1 
010 
10 0 


let the performance index to be minimized to be 


m 

l 


P 

I 


i=l j=i 


2 _ 2 
“i 


Solution; 

tt T - [9.33 13.33 -5.33] 



125 


and the modal controller feedback matrix is 


F « 

0 


4.66 

—4 .66 

• 6.66 

- 6.66 

2.66 

2.66 


6.6.3 Example 3 : Problem Considered by Sridhar and 

Lindorff [87] and Sirisena and Choi [88] 


A *» 


10 0 0 

0 2 0 0 

0 0-30 

0 0 0 -4 

]D 


B 


1 

0 

1 

1 


0 

1 

0 

1 


110 0 
0 0 11 


It is desired to shift eigenvalues (1= 1,2) of 
the uncontrolled system to locations ( p= -1, -2) without 

disturbing the eigenvalues ( 1 = -3> -4). 


Solution: 

In this case: 

p » 2, it* 2, m a? 2, r - 2, p + it _r = 2 > 0, 
kj** 6 and k 2 * -12. Taking c 2 = 1.0, we obtain from 

eqns. ( 6 * 18 ) 

« 1.0 

and •* [ l. 0 ] 
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Then gc T = [ 6 -12 ] 

and "the feed.bo.ck matrix is 

6 0 

-12 0 

The closed loop system matrix then is 

7 6 0 0 

-12 10 0 0 

6 6 -3 0 

—6 —6 0 — ‘ 4 

If now an attempt is made to shift the eigenvalues 
( A - -3» -4 ) to some other locations then it is found that 
the rank [D] = 4 and p+ H -r jf 0, which violates eqns. (6.19), 
These modes, therefore, now cannot be shifted by this 
technique. 

6,7 ACCESSIBLE STATE FEEDBACK DYNAMICAL MODAL CONTROLLER 
An interesting scheme for designing feedback 
controllers for systems with inaccessible state variables 
has been given by LioU et al [89]. The procedure enhances 
the dimension of the state space by considering the control 
inputs also to be state variables. In the augmented system 
thus obtained, the control variables are the derivatives 
of the original control variables. An advantage of this 
approach is that the original control variables will not 


[A + B jx |3^ D] « 
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have discontinuities while the system is in operation and 
th.e controller will he 'smooth acting' . This may imply 
longer controller life, liou et al [89] use infinite time 
optimal linear regulator theory for the determination of the 
feedback controller for the augmented system. This approach, 
however, may not he computationally attractive for higher 
order system [Section 3.1]. The formulation of Liou et al 
is used here, hut the solution to the problem is obtained 
by the application of modal control technique developed 
earlier. 


6.7.1 Development : 

The system (6.1) is augmented such that 

» 

z = A z + B w 


where 


z 

A 

B 


r T T-.T'. 
[ x u ]j 

~A B 

0 0 


f 


f 




T 


= mxm identity matrix and 


( 6 . 21 ) 


A and B are (n + m) x (n«) and (n«) .x » ■ system *md control 
matrices respectively of the augmented system. A modal 
control lav to achieve the desired closed loop poles may 
bo obtained by the procedures given in Sections 3.4 and 4.4 as 
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Eqn. (6.22) can "be rewritten as 
w = F x + H u 


(6.23) 


wliere 


F = [ F : H ], 

F s= (mxn) state feedback matrix and 
H = (mxm) matrix associated with, the feedback of u 
The equation (6.23) can also be expressed as 
u -Hu = + F 2 -2 


where 


e 2 ] . 


t 

r— ' vector of inaccessible states and 
x = (n-r) vector of accessible states 


£•1 


The'system eqns. (6.1) can be written as 


— '• — 

»1 


A 11 

j— 

C\J 

<4 H 


^1 


i — 
w H 

1 

ia 

— . — 


A 21 

A 22 


1 

ro 

i 

+ 

] 

B 2 1 


u 


0-21 A 21 


If the (n-r) x r matrix A a has rank r, then 
)- 1 exists. Eqns. (6.25) then yields 

Sl = a (x 2 - A 22 S 2 - b 2 a) 


where 


a 


i—l 


T 


= [ A 21 a 2 i p A 21 

, i r»A-ntroller law obtained from 
Then, the dynamic modal co 

eqn. (6.24) &&& ( 6 * 26 ) is 
( H 


U 


F 1 aB 2 ) u = F i a -2 " ^ P l a ^22 


- f 2 ] % 


(6.24) 


(6.25) 


( 6 . 26 ) 


(6.27) 
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This can. be rewritten as 

a + ®i i = » 2 k 2 + s 2 

where 


»1 = H - F 1 aB 2 ] , 

* 

G-2 = ^i a and 

& 3 = - C * 1^22 “ F 2 ^ 

6.8 EXAMPLE: DYNAMICAL MODAL CONTROLLER POR A 

POWER SYSTEM*; 


(6.28) 

(6.29) 

( 6 . 30 ) 


The power system considered in Section 6.3.3 is 
also considered here. The state variables A 1 ? ^ and A h 
are considered to be Inaccessible. To apply the b 
procedure of Section 6.7, the estate variables are 
rearranged in the state vector such that : 


T 

x ± = [Al* f h ] 

& 2 = [ n, Av f> Av g , Agj A gf? A6 ] T 


(6.31) 


The augmented system (6.21) has the following 


seven dominant eigenvalues 


0.0, 0.0, -0.014 + jo.7986, -0.0572, __0.0772 ± 30.1146 . 

A modal controller is to be designed to shift the above 
eigenvalues to the following locations 

-0.4 ± 30.915* -0.5, -0.6, -0.96 + j0.72, -1 . 

The set of dominant eigenvalues contains a pair of 
repetitive eigenvalues at the origin. The modal control 
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procedure of Chapter 3 cannot therefore be directly 
applied. To overcome this difficulty, either of the 
following approaches may be used: 


App ro ach I 

Observing the system and the control matrices of 
the augmented system, it is found that the following control 
law shifts the two zero eigenvalues to the locations -1.0 
and -0.5 respectively, without disturbing the rest of 
the eigenvalues: 


w 

where 



0 0 0 0 0 
0 0 0 0 0 


0 

0 


0 -1 

0 0 



The closed loop system matrix A 1 with the above 
feedback is 

= A + B Fq 

Then a control law i' 1 to shift the remaining 
dominant eigenvalues to the desired closed loop locations 
is found. The resultant modal control matrix F is obtained 
by summing the two feedbacks laws Fq and F^. It 
found to be 

273 2 1.2668 - 2.9498 - 10.7221 1.13 - 0.3545 1.3796 

3040 0.4029 - 0.8079 - 4.7426 - 1.15 4.0930 0 

0.6008 - 2.8396 0.5837 

- 8.0436 - 0.4792 - 1.7460 
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The dynamic modal control law utilizing feedback from 
accessible states is then given by 

- + G 1 £ = & 2 Sfc + & 3 2* 

where 

2.8396 -0.5837 

J0.4792 1.7460_ 

0.4552 -13.2406 0.0 0.0 0.0 0.0~ 

-6.0936 - 4.4660 0.0 0.0 0.0 0.0 
5.2889 -0.5472 -21.0894 -11.1235 1.1300 1.6776 
-0.5148 -0.2089 -6.9263 0.6320 -1.1522 0.5116_ 

Approach II 

The dominant eigenvalues are divided into two 
groups, such that each group has only one zero eigenvalue: 

Group I : [ 0.0, -0.0114 + j0.7986, -0.0572 ] 

Group II : [ 0.0, -0.0772 + 3 0.1146 ] 

A modal controller feedback matrix P- L is obtained 
to shift the eigenvalues of group I to the desired locations. 
Then a feedback matrix 3? 2 computed to shift the eigenvalues 
of group II to the assigned locations. The resultant feed- 
back matrix F is obtained by adding the two matrices 
and P 2 . The matrices and $ 2 are selected to minimize 

the performance index 

m n 9 



s= 

? 
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where f ± . is the element of the feedback matrix F 
corresponding to i un row and j xxl column. 


The result and modal control law for the augmented system 
is found to be 


w = 


2.2407 0.8016 -2,9995 -1.5569 0.2332 -0.0700 

-3.8417 0,1574 -0.6008 2.7968 -2.0503 4.3960 

0.5978 0.5219 -3.2861 -0.0041 

0.6940 -8.9629 4.6945 1.2995 


z 


and the dynamic modal control law utilizing feedback from 
accessible states is 


u + G-^ u 


G 0 X-i + G~r X- 
2 -g 3 -^2 


where 


G-i 


Gr> 


Gi 


3.2860 0.0041 

-4.694 5 1.2995__ 

0.4552 -13.2406 0.0 0.0 0.0 0.0 

—6 . 0936 —4.466 0.0 0.0 0,0 0.0 j 

5.2889 -0.5472 -21.0894 -11.1235 1.1300 1.6776 

-0.5148 -0.2089 -6.9263 0.6320 -1.1522 0.5116 

The closed loop system responses with the above 


controllers with the initial conditions considered in 
Section 3.6.i are given in Figure 6.3 . Ihe improvement 
in system performance is clear from these figures. 



With modal controller 
Sec» 6.8 .Approach I 

With modal controller 
Sec. 6.8 Approach II 


Tome in Sec. 
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6.9 CONCLUSION 

Techniques for designing accessible state feedback 
and output feedback modal controllers for linear, time 
invariant, multi-input-multi-output systems have been 
developed. Procedures for optimization of controller 
gains have been indicated. The assumptions made in the 
development of the algorithm do not seem to be overly 
restrictive. A method for designing dynamic modal controllers, 
based on the work of Liou et al [89], has also been developed. 
This is computationally more attractive than the procedures 
given in reference [89], especially for high order systems. 
Numerical examples illustrating the application of above 
procedures have been given. 



CHAPTER 7 
CONCLUSIONS 

7.1 GENERAL 

The basic aim of this thesis has been to develop 
suitable techniques of designing controllers which improve 
the dynamical behaviour of power systems. These 
techniques should be such that the designer has a 'feel* 
for the problem and his engineering judgement can be 
better utilized. The work reported and the contribution 
made in this thesis are reviewed in this chapter. Suggest- 
ions and scope for future work in this area have also been 
indicated. 

7.2 REVIEW OF THE WORK DONE 

Following a brief survey ox the techniques available 
in the literature for designing control schemes for improv- 
ing the dynamic performance of power systems, the state 
space descriptions of the power systems under invest igation 
are discussed in Chapter 2. In Chapter 3, a procedure for 
designing regulators based on modal control theory with 
an optimality criterion has been suggested to obtain 
improved performance of linear time invariant dynamical 
systems. The algorithm developed is completely general, 
independent of the system structure and computationally 
easy. However, the dyadic structure assumed for modal 
controller matrix may be restrictive and may result in 
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large feedback gains especially in high order systems 
with large number of inputs. An algorithm based on the 
grouping of the dominant eigenvalues is therefore proposed 
in Chapter 4 resulting in lower feedback gains even for 
high order systems. A criterion for dividing eigenvalues 
into groups has been given. Procedures for the design of 
sub-optimal and near— optimal modal controllers requiring 
smaller computational efforts have also been developed. A 
technique is suggested to identify dominant control inputs 
for assigning the eigenvalues in a particular group. The 
feedback gains are marginally affected when feedback is 
arranged to dominant inputs. The techniques developed in 
the chapter are demonstrated for a large dynamical system 
described by fortyone differential equations. An effective 
method based on eigenvalue sensitivity to design modal 
controllers for improving the system performance over a 
wide range of operating conditions is developed in Chapter 
5 . Thjfcs wide range modal controller reduces the variations 
in the closed loop system eigenvalue locations with the 
change in operating condition. Techniques for designing 
accessible state feedback and output feedback modal controll 
ers for linear time invariant multi-input multi-output 
systems have been presented in Chapter 6 . Procedure for 
optimization of controller feedback gains is suggested. 

A method for designing dynamic modal controller is also 

reported. 
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7.3 SUGGESTION AND SCOPE FOR FUTURE WORE 
7 « 3*1 Assignment of Zeros in Addition to Poles 

In this thesis the improvement in the power system 
performance is obtained hy assigning the dominant eigen- 
values to favourable locations. However, there is no 
control on the location of zeros of the closed loop system 
which also affect the system dynamic performance. Therefore 
it is desirable to develop techniques to assign both the 
poles and zeros to favourable locations. 


7.3.2 Mult ima chine Sy st ems 

Examples of single machine connected to infinite 

bus are considered to illustrate the techniques developed 
in this thesis. It will be worthwhile to consider a 
realistic multi-machine power system to investigate the 
suitability of modal control techniques. 


7.3.3 Attainable Pole Placement with Output Feedback 

The algorithm developed to design output feedback 
modal controller for linear systems requires certain 


conditions to be satisfied. It may not always be possible 
to meet these requirements. Therefore it is desirable to 
investigate the attainable pole placement positions with 


a given output matrix. 


7 . 3*4 


Decentralized Control 
The modal control law requires 


the feedback of 
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all the variables* In a modem power system which consists 
of large number of interconnected subsystems (Areas) it may 
be difficult and uneconomical to arrange feedback from 
the other areas to a controller in a particular area* It 
is worthwhile to investigate whether it is possible to 
improve the system performance by local control, that is, 
information available in a particular area is fed to its 
control inputs. 

7*3*5 Non-linear Model 

The performance equations are linearized around an 
operating point and control scheme is designed to improve 
the performance of this system for small disturbance. It 
is desirable to investigate the performance of non-linear 
model with the linear regulator for large disturbances . 

7*3.6 Stochastic Control 

In the studies reported in this thesis, a disturbance 
of deterministic type is assumed. In practice the type of 
disturbances are usually stochastic in nature. The problem 
of stochastic control with state dependent noise and 
unknown noise statistics is worth studying* 
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APPENDIX A 


MODAL control 
C onsider the system 

i = A x + B u (A.l) 

and its related single input system 

x = A 3C + a z (A. 2) 

where vector a is generated from the columns of the control 
matrix B such that 


m 

a = l a. b. = B a (A. 3) 

- j =1 3 “3 


where d.-'s are scalars and 
control matrix B. 


h.'s are the columns of the 
3 


Consider the transformation 


x 



(A. 4) 


where ¥ = [w^ w £ f ...» wj, the eigenvector matrix A. 

Let V = [vj_, v 2 , ...» Le the reciprocal eigenvector 

matrix of A and the eigenvectors are normalized such that 

V T = W” 1 U-5) 

Then 

A = V T A W 

- Dia[ X 2 » •"» (A.6) 

Using the transformation (A.4), the system (A.2) is trans- 
formed as 
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y = A£ + V T a z 


Ay + p Z 


(A. 7) 


m 

where p = v a 


Consider a feedback control law of the form 


z = 


a 

i 

i=i 


h y i 


l 

l k < v, , x > 
i=l 1 ^ L 


T 

g £ 


(A. 8) 


where the gain vector g = I k. v. and k *s are the 

i=l 

controller gains and v\ 's are used as the measurement 
vectors. 

The closed loop system is 


L 

y =/ A'y 


where A = ,A + £ k 


T 


k = {k-,* k , ...» k , 0, 0] T is ann-vector 


1 2 

and k ± is the controller gain and 


A 


A 11 -° 


A A 

21 22 
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where d. - = p. k. 

13 i o 

if d . . = 0 , then the if* 1 mode is uncontrollable. 

-** J 

The characteristic equation Det! si - A i = 0 decomposes into 

_ n 

Det | si - An I n (s - X. ) =0 (A. 10) 

i=4 +1 


Thus the first eigenvalues are variant and the rest 
(n-£) eigenvalues remain invariant. If the eigenvalues 
of the controlled system are (p-^, P ^ * • • »P^ »^^ + i» • * * »A n )» 
then 


Det f si - A i;L I = (s-A^Cs-Pg) ... (s-P £ ) = 0 

$et j A-I- A 2i | = ( *0 = * f 3 ~ l*2» • • * »*■ ) (-4.11 ) 


or J | ~~ j P j 9 ^ if 2 $ • * • fXt } (j 

The determinantal equations after simplification yield 


k = [ ik (P . ] /CP-? ® (3""1> 2f • • • jfc ) (A#12) 

j L i 3 lid 1 3 

±¥i 
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Substituting p.. = v^ £ and a = B a, the eqn: (A.12) reduces 
to 

a % 

k 3 = c Ji ( p i- V ]/[s j B - •! < W ] U.13) 

Vi 

and 

k d + l = k ;j lf *5+1 = A 

The gain vector g and the control law z are realized from 
eqn.(A.8). A real control law is obtained for shifting the 
complex eigenvalues in conjugate pairs. The control law for 
multi-input case is then realized as 

u = a z (A. 14 ) 

and 

u ^ F x 

Then the modal controller matrix 


T 

F = a g 


(A.15) 



APPENDIX B 


LAGRANGE MULTIPLIER TECHNIQUE 


UasQ l„i . Let the function J to be minimized is Jg (eqn.3«10)* 
m n 

J = I a. + E ff? C'g i \ 

3=1 3 1=1 1 l J 

and the equality constraints to be satisfied are 


g 


H 

l k t Zi 
i=l 1 1 


0 


where v^ and are defined in (A.13)* 


(B.2) 


Using Lagrange multiplier technique for the optimi- 
zation problem, the Lagrangian function is formed by 
augmenting the function J with the equality constraints 
multiplier by the respective Lagrangian multiplier 


that is 

= 


m n m & 

•p = i «? + I gf + HLg - 1 k n X-;] 
1 j=l 3 i=l 1 ~ ~ j=l J J 


(B.3) 


\ + l g| + I^Le - 
3=1 3 i=l 1 D= 

where the elements |i^» (i = l»2,«*.»n) of the ]i are the 

«n' Lagrange multipliers. It follows from eqn.(B.5) that 

the optimum value of £ f , when the equality constraints 

(eqn. B.2) are satisfied, is equal to the optimum value of 

J. Then it follows 


Tg7 


<3 


SS fg. — £ k - *4 

1 j=l 3 3 


o 


(i * 1,2,...,: 


(i = 1, 2, • . • , 


(B.4) 

(B.5) 



152 



where 


2 “ 1 + 11 [ jix [( P ="i 

C ( n (V x i»(4 E ^^> 

p=l r j j 

(i = 1,2, ... ,m) (B»6) 


is the element of v^ and 

is the i^^ 1 column of the control matrix B • 


At optimum value of J and£ f e qns. (B.4) to (B.6) must be 
equal to zero. The negative gradient of the Xagrangian 
function£. f with respect to a is 

( A a = - (3ty3a) (B.7) 


Case 2 : let the function J to be minimized in J 1 (eqn.3«8). 


m n 

J = Z E ?• g- 
3=1 i=l 3 


Then the lagrangian function is 



m n 

Z S 
0=1 i=l 


a? 


4 * j?[« - 4 k 3 2 s ] 


Then it follows 


3 ^f 

3^ ± 


= ig* - 


a 

z 

d=i 


k. v? 1 ^], (i = lj2,.*»»n) 
3 3, 


(B.8) 


(B*9) 


(B.IO) 



-p m . 

= E 2 a “ g + (i, , (1=1,2, ...,n) (B.ll) 

5 Si 3=1 3 1 1 

E 2a ± g 2 + / [ * [( " ( p p a 3 ))p 3i S ]/ 

3=1 J D=1 P— L 

[( n (>■ -^))(^ b a) 2 33 

P =i p 3 J 

(i = 1,2, (B.12) 

The negative gradient of the lagrangian function f f with 
respect to a is 

( A a) = -(<* %/^0 * 



(B.13) 



APPENDIX C 


LARGE OEMICAL PLANT 


A generalised chemical plant model of thirty first 
order was postulated by Williams and Otto [78] for the 
investigation of computer control. Davison and Chadha [73] 
modified the above model by increasing its order to forty 
one by considering four plates of the distillation col umn s 
instead of two plates considered earlier to give a realistic 


model. The description of the plant as well as that of 
linearized mathematical model with data is given in detail 
in reference [73]# 

BRIEF DESCRIPTION OF PLANT 

The basic units of the chemical plant are a reactor, 

a heat exchanger, a decanter and a distillation unit complete 
with reboiler and condenser. Two pure feeds and one recycle 
feed are provided to the reactor. Six materials are present 
in the effluent from the reactor. The exothermic reactions 
are controlled by the addition or substraction of heat 
to the reactor. 


Each of the reactor and heat exchanger units are 
represented by a set of six differential equations in the 
mathematical model. five differential equations are required 


to represent 


unit and the distillation plant Is 
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represented by a set of twenty differential equations. 

The interconnections are modelled by four differential 
equations. In all forty one differential equations represent 
the chemical plant. 


DATA FOR THE LINEAR MODEL OF THE PLANT 
The data for the linear model 
x = A X' ■+. B u 

as given in reference [73] are as follows: 

The system matrix A is given by: 

A* ( a^ ,. ) , i=l ,2, • • • , 41 

3 = 1 , 2 , . » « , 41 

a l,l “ ~4°*74, a i f 2 ~ a- 2f i = ”20.1, a 2t 2 ~ 

a 2 3 ” a 3,l = ^°* 2 * a 3,2 = ^*^ 8 * a 3,3 = *“ 182, 44» 

a 


a, 


3>5 = -7.35, a 4j2 = 8.94, a 4>3 = 136.8, a 4>4 = -20.64, 

3>2 = 4.47, a 5f3 = 55.9, a^ = -24.315, a gj3 = 37.5, 
a g 5 = 11.025, a g>6 = -20.64, a ? ^ = 3030.0, a^ 2 = 1449.5, 

= 436.0, 


= 8168.0, a ?>5 = 190.6, a ?>7 = -35.82, a 8>? 
to o = -603.4, a Q>9 = 167.0, a g>8 = 17.7, a g>g 


= -202.69, 


a 


21,8 


= 16.42, a 21,21 ~ " l6 '^ 2 * a 22,22 


= -1.56, 


a 23,23 = - 1 - 04 ’ a 24,24 
a 22}2? = 1 - 56 ’ *23,28 = 1 * 04 ' *24,29 


= -2.08, a 25j 25 = -3* 6 4, a 26,26 — 3 ‘ 64 * 

25,30 


= 2.08, in — 3.64 j 


'26,31 


3.64 



and the following algorithms give the rest of the 
non-zero elements. 


a i-36,i = 10 -75 

f 

1=37,38, .. ,41 

a i,i “ -350.0 

f 

1=10,11,.., 15 

a i+9,i = 350.0 

t 

1=1, 2,.., 6 

a i+6,i = 1 3' 7 

9 

1=10, 11,.., 14 

— 



. = —13 #2 

J-fX 




f 

i = l6 jlT f « • f 20 

a i+ll , 1 = 15,6 



a i+5,i+5 = “ 46, ° 




a i +10,1+10 = “ 55,95 

a i+15,i+15 = ~ 55,65 


a i+5,i = 


14.8 


a i+5,i+10 = 15,6 
a i+10,i+5 = 55,65 
a i+10,i+15 = 18f5 
a i+15,i+15 = 55,65 


1=22, 23 ,.., 26 


1 = 22 , 23, *• *26 


The control matrix B is given by 
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r5 


and its non-zero elements are: 
b l,l = 2.16x10 4 , b 2j2 = b Xjl , b lj8 = 2.87 x 10 
b ?>1 = 1.15X10- 2 , b ?>2 = b ?A , b?j3 = 1.5, b ?>4 = 1.5, 

b 7 >5 = 3,36 » b 9,6 = 185 * 0 > b 9>? = -2.4x10-2, b 2>8 =9.07xlO- 5 , 
b 3,8 = 0-582xlO“ 5 , b 4>8 = 8.21 x 10“ 5 , b 5j8 = 0.819 x 10“ 5 , 
b 7,8 = ~1 345.0 x IQ" 3 . 



APPENDIX T) 


EIGENVALUE SENSITIVITY [80] 

let f(s,A) be the characteristic polynomial of A. Then 
f ( S » A 0 ) = det [if - A 0 ], 

Consider the variation dA in the nominal value A an d let 

o 

the eigenvalue differential be ds, due to differential change 
dA. 

On expanding f(s,A) = 0 in Tayler series about nominal 
value A , we obtain ©© 

Q /-v 

f (s» A) * [f(s,A)] 0 + l [d k f(s,A)] 0 = 0 (D.l) 

k— i 

k fh 

where d (.) is the k differential of (.)* 

Since [f(s,A)] Q = f(s,A Q ) = 0, the eqn. (D.l) reduces to 

f(s,A) = l [d k f(s,A)] 0 = 0 (D.2) 

k=l t 

Erom matrix theory [90'], it can be shown that for any square 
matrix P 

d [det P] = adj P 9 dP (0.3) 

sud the inner product of two e^uidimensionaX square matrices 

is defined as 

e» , T r i v. ^ V 

0 | D a Z c • d. > v w 
ill 

where c. is the i th row of 0 and d ± is the i th column of B* 
i 

Then from eqn. (B.2) 

^ [d^ 1 [adj (sI-A) ^d (Is-A)]] 0 * 0 


(D*4) 
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If the differential change is linear, then a k A - 0 for 

1 and if the variation in parameters are sufficiently 
small *th.en the terms 


oc 

k =2 “ET C4 k_1 [adj(sl-A) « d(Is-A)]L 

are small being the product of differentials. Then the 
6C|ii* (D*4) is rewritten as 

[adj (sI-A) ] Q « Ids = [adj (sI-A)] 0 © dA ( 

Since adj (sI-A) © Ids = Tr.[adj (sI-A)J ds, where Tr.( . ) 
is the trace of the matrix (.), the eqn. (D.5) reduces to 


Tr. [adj (sI-A)] q ds = [adj (sI-A)] q © dA 
and for a particular s^. the above equation yields 

d®i = [ ad 0 ( sX— A) J q J ■*" [[adj (sI-A)] o © dA] 

= [Ir. E°(A i d )]- 1 [E° (X^)^ # [a (A + b f, 1 )] 

ds. - —1 d m 

or d©” = ^ Tr * R °^ x i^o x i^o # ( A + & f 0 )L 

j ^''j - u . 

( i~l # ♦! j n) f ( j —l f 2 f % • • ,r) 

where R c ( \*) = adj [ X* I - A ] and 
d th 

A ^ is the i desired closed loop system eigenvalue. 


(D.6) 
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